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Abstract 

The concept that a chemical reaction may go in the direction of increased total equalized 
electronegativity (Xeq) has been explained for substitution reactions involving complexation by 
monodentate ligands and chelation by polydentate ligands. It is observed that increase in Aeq is even 
more in case of substitution reactions involving polydentate ligands which leads to products having 
greater number of chelate rings. 

(Keywords : electronegativity/chelation/stability of metal chelates) 


Introduction 

Although the principle of electronegativity equalization 1 has been cited in chemistry 
literature since 195 1 its application [equalized electronegativity (Xeq)] has not been fully 
exploited except to calculate the partial charges on atoms and groups 2 , and to Pearson 
hard-soft-acid-base reactions 3 . Chelating ligands, the chelate effect and the stability 
increment resulting from chelation is well known 4 . The chelate effect for a n-dentate 
chelating ligand refers to the reaction where the chelating ligand replaces the n-unidentate 
analogous. 


M (unidentate) n + Chelate ^ M (Chelate) + n-unidentate 
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In this work, some applications of Xeq to substitution reactions involving complexation by 
unidentate ligand and chelation by polydentate ligands are discussed. 

Method 

Xeq values for various species in substitution reactions were calculated using Bratsch 
method 5 according to equation 


Xeq = 




where .Y = v = the number of atoms in the species formula, q is the charge of the species 
and X is the Pauling electronegativity. 

Results and Discussion 

The following examples (Table 1, reactions 1*13) of the substitution reactions are now 
examined in terms of total Xeq values of the products and reactants. Reactions 1-13 are 
substitution reactions of the type 

M (H 2 0)+ m + L' n = M (H 2 0) x ^ L +(m ' n) + H 2 0 

If L were a polydentate ligand then two or more molecules of water would be displaced 
from the coodination sphere of the complex. It is obvious that Xeq i.e. (2Xeq Products - 

.wA'eq Reactants) is greater than zero for reactions (1- 13), thereby suggesting that a 
substitution reaction of above type proceeds in the direction of increase in the total Xeq 
\ time of products in comparison to the total Xeq of the reactants. Some results are compiled 
from literature on stability constants of complexation reactions 1-13 (Table 1). It is observed 
that increase in Xeq value i.e. AXeq is greater for reactions 2, 4 and 6 involving bidentate 
rigaud <en) as compared to AXeq for reactions 1, 3 and 5 respectively involving only 
unidentate ligand. The increase in AXeq between nonchelated and chelated species forming 
reaction* is 4.79 (reactions 1 and 2) and 2.31 (reactions 3 and 4; reactions 5 and 6). This 
increase in AXeq may be related to the stability of the chelated species. 

It is evident from log values (reactions 4, 7 and 8, Table 1) that stability increase 
with more substitution of H 2 0 by en from coordination sphere of the metal ion. These 
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Tabi 

le 1 - AAeq values and stabiht} data for substitution reactions in< 
species the number refers to its Xtq value. 

*oiv ;ng compieud 

on*'che 

labor., Below each 

Sr. 

No. 

Reaction 

Aec 
' Z.Vec 
Products- 

V y A( , 

Reactants 1 

i A 

* log 0, 

ly Constant 

I. 

Cd(H 2 0)T + 4CH 3 NH 2 = [CdtCH 3 NH 2 ; 4 f + f 4H ; 0 

2.78 2.34 2.36 2.52 

0.22 

log 

= 6.5 (Ref. 7) 

2. 

CdiH 2 0)T + 2 en = [CJ.er. : ; 2 ' + 4H-0 

2.78 2.36 2.51 2.50 

5.01 

log fin 

= 12.6 iRef. ?j 

3 

-t- 2NH 3 = ;NuH : 0. 4 ;.NH 3 ;2: 2t + 2H : 0 

2 72 2.36 2.64 2.50 

0. 2d 

»Cg 32 

= 5.0 . Ref. S,» 

4. 

Xi(H 2 0}^ h- en = [NuerOiH 2 Qi 4 r + t 2H : 0 

2.72 2.36 2.59 2.53 

2 51 

log 2; 

= 75 Ref. 8 

5. 

Cu.H 2 0,:* + 2NH 3 = LCu.H : 0 4 NK 3: ; 2 ‘ - 2H ; 0 

2.73 2.36 2.65 2.50 

0.23 

log p 2 

= 7.7 .Ref. 9* 

b 

Cu\H 2 0;£ T + en = ;Cu(H 2 0i 4 en’.' + •* 4H 2 0 

2 73 2.36 2.60 2.50 

2.51 

log 3: 

10.6 i Ref. 9) 

7 

N:.H 2 0;5* + 2er. = [N\H 2 0 i 2 ien) ; T + 4H ; 0 

2.72 2.36 2.52 2.50 

5.08 

iog p 2 

= 13.50 (Ref. SO) 

S 

Ni(H : 0)g + + 3 en = [Ni(en) 3 ] 2+ + 6H ; 0 

2.72 2.36 2.47 2.50 

7.67 

’•og P 3 

= 17 6 (Ref. 10) 

9. 

Ni(H 2 0)g + i- Gly = (Ni(H 2 0) 4 Gl> ] + + 2H 2 0 

2.72 2.28 2.61 2.50 

2.61 

log Pi 

= 6.2 (Ref. S; 

10. 

Ni(H 2 0)g + + IMDA 2 " = [Ni(H 2 0) 4 UMDA)] + 3H 2 0 

2.72 2.25 2.54 2.50 

5.07 

log P| 

= S.O (Ref. 8) 

11. 

Ni(H 2 0) 2+ + NTA 3 " = [Ni(H,0) 2 (NTA)f + 4H 2 0 

2.72 2.24 2.48 2.50 

7.52 

log Pi 

= 1 1.3 (Ref. 3) 

12. 

Ni(H 2 0)g + + dien = [Ni(H 2 0) 3 (dien)] 2+ + 3H 2 0 

2.72 2.36 2.54 2.50 

4.96 

log Pi 

= 10.5 (Ref. 8) 

13. 

Ni(H 2 0)5 + + trien = [Ni(H 2 0) 2 (tnen)] 2 + + 4H 2 0 

2 72 2.36 2.50 2.50 

7.42 

log P[ 

13.7 (Ref. 8) 


^Abbreviations : en-H 2 N CH 2 CH 2 NH 2 ; pn - H 2 N(CH 2 )3 NH 2 ; dien - H 2 N CH 2 CH 2 NH CH 2 CH 2 NH 2 ; 
irien-H;N CH 2 CH 2 NHCH 2 CH 2 NH CH 2 CH 2 NHy, IMDA 2 ' - HN (CH : COO'h; NTA 3- - N'(CH 2 COO') 3 
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observations find support from AYeq values which increase on substitution of H,0 by 
greater number of en. 

It is well established that a ligand with greater number of chelate rings forms even the 
more stable complex 6 . This is obvious from stability data (log P, values) for a series of 
amino-carboxylate ligands GLY' < IMDA 2 " < NT A 3- (reactions 9, 10 and 1 1 respectively). 
Likewise in the amine-ligand series, the stability constants increase as : en < dien < trien 
which correspond to reactions 4, 12 and 13 respectively. Such stability order is evident 
from AXeq values (Table 1) which also follow the sequence : GLY - < IMDA 2 ' < NT A 3- ; 
en < dien < trien. The number of chelate ring formed is one each in gly" and en, two each 
in IMDA 2 ' and dien, and three each in NTA' 3 and trien. It is, therefore, evident from AZeq 
values that more is the number of chelate rings formed by a poly dentate ligand, the larger 
is the increase in Meq in substitution reactions involving chelation. 
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Abstract 

New synthetic azo reactive dyes have been prepared by coupling diazotised 4,4’ -methylene 
bis-orthoanisidine with various anilino cyanurated coupling component and their dyeing performance 
as hot brand reactive dyes have been assessed on viscose rayon, silk and wool fibres. These dyes have 
been found to give a wide range of violet to brown shades with very good depth and levelness on 
each fibre. These dyes were characterised by elemental analysis, IR and PMR data. 

(Keywords : IR/4,4’-methylene-bis ortho anisidine/bis-azo reactive dyes/viscose rayon/ 
silk/wool). 


Introduction 

The mono chloro triazinyl reactive system is of major commercial importance in reactive 
dyeing. 1,2 In recent years a large number of monochloro triazinyl reactive dyes of different 
types have been synthesised for natural and synthetic fibres 3,4 . The research based on 
bifunctional reactive dyes has been widely considered due to there higher fixation yield on 
various fibres 5 . The dyes with two reactive groups give a higher fixation than dies with 
one reactive group 6,7 . In the present paper, we report the synthesis and dyeing performance 
of ten hot brand bifunctional reactive dyes. 

Materials and Method 

All the melting points (m.p.) were determined in open capillaries and are uncorrected. 
The IR spectra were recorded on Shimadzu Japan DR 8001 (FTIR) spectrophotometer by 
using Nuzolmull techniques and PMR spectra on a Varian 300 MHz instrument using TMS 
as internal standard and CDC1 3 + DMSO-d 6 as solvent. Fastness tests were assessed by the 
standard method of testing (BS : 1006 - 1978 and IS : 765 - 197). 
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Preparation of 4,4' -methylene-bis-ortho-anisidine (A) : Ortho-anisidine (12.3g, 0.1 mol) 
was dissolved in water (125 ml) and 36.5% hydrochloric acid (25 ml) at 50°C. The reaction 
mixture was then reacted with 3% aqueous formaldehyde (25 ml) solution at 60°C with 
stirring for an hour and neutralized with 10% sodium hydroxide, when yellow precipitates 
were obtained and filtered, washed with hot water, dried and crystallised from acetic acid, 
Yield 88%, m.p. 100 - 101’C (Found : C, 69.40%; N, 10.60%. C 15 H 18 0 2 N 2 requires : C, 
69.76%; N, 10.85%). 

Teirazotization of 4 , 4 ' -methylene-bis-ortho-anisidine (B) : 4,4’ -Methylene-bis-ortho- 
anisidine (A) (2.58 g, 0.01 mol) was suspended in H 2 0 (60 ml). Hydrochloric acid (10 ml) 
was added dropwise to this well stirred suspension. The mixture was gradually heated upto 
?0*C. till a clear solution was obtained. The solution was then cooled to 0-5°C in an ice-bath 
and a solution of NaN0 2 (1.2 g) in H 2 0 (8 ml) previously cooled to 0°C, was then added 
to it over a period of 5 min. with stirring. The stirring was continued for an additional hour, 
maintained at the same temperature, with positive test for nitrous acid on starch-iodide 
paper. After destroying excess of nitrous acid with required amount of a solution of 
sulphamic acid, a clear tetrazo solution was obtained which was used for subsequent 
coupling reaction (Reaction 1). 

Prepartion of anilino cyanurated coupling components (R) : The coupling components 
were prepared in two steps (Reaction 2). 

(i) Preparation of cyanurated coupling components (R") : A solution of cyanuric 
chloride (3.70 g, 0.02 mol) in acetone (50 ml) was stirred at a temperature below 5°C for 
a period of an hour. A neutral solution of coupling components (R’) (0.02 mol) in aqueous 
sodium carbonate solution (10% w/v) was then added to it in small lots in about an hour. 
The pH of the solution was maintained at 7 by simultaneous addition of sodium carbonate 
solution (1% w/v). The temperature was maintained below 5°C for further 4 h when a clear 
solution was obtained. The cyanurated solution (R") thus formed was used subsequently. 

( ii) Formation of anilino cyanurated coupling component (R) : The cyanurated coupling 
component (/?") was heated upto 40- 50°C for half an hour. Then aniline (1.86% g. 0.02 
mol) was added to it dropwise at the same temperature, during a period of 30 min. 
maintaining the pH at 7 by simultaneous addition of sodium bicarbonate solution (1% w/v). 
After the addtion was completed, stirring was continued for further 3 h. The anilino 
cyanurated coupling component (R) thus obtained was used for further coupling reaction. 

Formation of Dyes (D , to D 10 ): To an ice-cold and well stirred solution of anilino cyanurated 
coupling component (/?), a freshly prepared solution of tetrazo soluton ( B ) was added 
dropwise overa penod of 10-15 min. The pH was maintained at 7.5-8.5 during the addition 
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(A) 


NaNo/HCI 


Diazotisation 


M/ 


h 3 co och 3 

/ r >-N=N-C! 


(B) 



Where, R = Various anilino cyanurated coupling components 
for the formation of D , to D 10 (Table - 1). 


Reaction 1 
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S0 3 H sq 3 h 



Reaction 2 



Table 1 - Yield, A, max , m.p., Rf and elemental analysis of the prepared dyes 
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and stirring was continued for 4 h, maintaining the temperature below 5°C. Sodium chloride 
(12 g) was then added and the mixture was stirred for an hour. The solid dye separated out 
was filtered, washed with minimum amount of acetone and dried at room temperature 

(Table 1). 


Results and Discussion 


The dyes were violet to brown colour and obtained in excellent yield (72 to 86%). The 
purify of all dyes has been checked by thin-layer chromatography 8 . The yield (%), m.p., 
/.max, Rf value and elemental analysis data of dyes are reported in Table 1. The shade and 
the fastness properties are reported in Table 2. 

IR Spectra : The IR spectra of dyes D, to D 10 showed characteristic band at 2940-2925 
(-CH,-) 9 , 1730-1725 (C=0), 775-760 (S-triazine), 725-715 (C-Cl), 1480-1470 (-NH), 
1385-1380 (-N=N-), 1090-1075 (S=0), 3450-3400 cm" 1 due to O-H stretching vibrations. 
The IR spectra of some hot brand reactive dyes are shown in Fig.l, 2. 

PMR Spectra : In general the dyes show a singlet peak at 8 4.80-4.05 (2H, -CH -), 
5.78-5.08 (2H, Ar-OH), 4.80-3.78 (2H, -NH-Ar) and around 7.5-8.7 (aromatic protons) 

Dyeing : The dyes varied in hue from light to dark shade of violet to brown colour 
depending on the nature of substituent present in the coupling component. The result shows 
that the presence of the e~ donating groups increase the bathochromic effect. However, the 
light fastness appears to be good to excellent because of higher molecular weight. The 
overall dyeing properties appear to be good to excellent. 
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Abstract 

The interaction of Mg(Il), Mn(II) and Cd(II) ions with cephalexin, a wide spectrum antibiotic drag 
of cephalosporin group, in presence of biologically important ligand a-alanme has been inv estigated 
by potentiometric studies. The proton dissociation constant of cephalexin, stability constants of the 
binary complexes of 1:1 and 1:2, and ternary complexes of 1:1:1 molar ratios of the metal ion to 
ligand have been reported at 25’ and 35’C and 0.1M (NaCl 4 ) ionic strength. The order of stability 
constants follows as Cd(II) > Mg(II) > Mn(ll). Thermodynamic parameters |AG, AH and AS) 
calculated at 35*C and 0.1M ionic strength indicate that the drug-metal ions and alanine interactions 
are enthalpy characterised. 

(Keywords : metal-drug equilibria/stability constant/mixed-ligand complexes/Cdi II )/Mg(Ii) 
complexes/cephalexin complexes) 

Introduction 

It is well known that ternary complexes play an important role in biological processes 
in which enzymes are known to be activated by metal ions 1 . Ternary complexes have also 
been implicated in the storage and transport of active substances through biological 
membranes. These phenomena depend upon the formation of higher order complexes and 
the nature of involed metal ions with antibiotic drugs. In the presence of amino acid the 
enzyme may yield useful information for elucidating the molecular mechanism of drug 
action. In this paper, the results of the study of metal-cephalexin-alanine equilibria in 
solution are reported. 

Cephalexin is one of the most important semi-synthetic cephalosporins. Cephalexin 
has benzalamine (C 6 H 5 CH.NH 2 ~) and hydrogen (H) groups present at R t and R 2 respectively 
in 7-amino-cephalosporanic acid molecule. 

Materials and Method 

All the chemicals used were of A.R. or S.M grade. Cephalexin sample was obtained 
from Lupin Pharma Ltd., Mandideep, Bhopal (India). 
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O 



H H H A 

-vi r c ' 


CH. 



COOH 


(I) 7-Aminocepha!osporanic acid 



(Jl) Cephalexin 


The preparation and standardisation of the metal ions and the ligand solutions, and the 
experimental details are the same as reported earlier 2 . 

Stoichiometry of the complexes formed was ascertained by conductometric titrations. 
The Calvin-Bjerrum’s pH-titration technique as modified by Irving and Rossotti 3 ' 5 was 
applied to determine various formation constants at 25° and 35°C and 0.1M (NaC10 4 ) ionic 
strength. For pH-titrations the following thermostated mixtures were titrated with the 
carbonate-free 0.1M NaOH solution: 

(a) 10 ml 0.04M HC10 4 + 1.0M NaC10 4 

(b) Mixture (a) + 5 ml 0.004M cephalexin 

(c) Mixture (a) + 5 ml 0.004M alanine 

(d) Mixture (b) + 10 ml 0.002M metal 

(e) Mixture (c) + 10 ml 0.002M metal 

(f) Mixture (b) + 5 ml 0.004M alanine 


(g) Mixture (f) + 10 ml 0.002M metal 
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Table 1 - Stability constants and thermodynamic parameters p = 0. 1 M (NaCC>4); Temp. = 25°C (a) and 35“C (b) 
x = half integral, y = Least square method. 


Constants 

H 

Mg 

Mn 

Cd 

logtff 

(a) 3.50 y 





(b) 3.25 y 




log /if 

(a) 8,75 y 





(b) 1.50* 




4 l 

(a) 

5.35 x 

4.47 x 

5.77 x 



5.27 y 

4.47 y 

5.60 y 


(b) 

5.00 x 

4.25 x 

5.22 x 



5.00 y 

4.21 y 

5 . 20 * 

*°g ^ ml 2 

(a) 

3.45 x 

3.37 x 

3.55 x 



3.40 y 

3.32 y 

3.55 y 


(b) 

3.32 x 

3. 15 x 

3.45 x 



' 3.2& 

- 312 y ... ■ 

■ , '-3.42 ifc: *‘-' ' 

iogp 2 

(a) 

8.70 x 

7.84 x 

9.32 x 



8.67 y 

7.79 y 

9.15 y 


(b) 

8.32 x 

7.40 x 

8.67 x 



%.2<? 

7.33 y 

8.62 y 


(a) 

4.93* 

4.38 x 

5.28 x 



4.96 y 

4.36 y 

5.25 y 


(b) 

4.58 x 

3.95 x 

4.77 x 



4.55 y 

3.92 y 

4.78 y 

Alog K 





lo 8 ~ *MlJ 

(b) 

-0.45 y 

-0.29 y 

-0.42 y 

-AC (Kcal/mole) 


6.46 y 

5.56 y 

6.78 y 

-AH { — do — ) 


17.31 y 

18.57 y 

19.84 y 

-AS (cal/deg/moie) 


+35.22? 

42.24 y 

42.40 y 
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The ionic strength of the above solutions was maintained to 0.1 M by adding required 
quantity of LOM NaC10 4 solution and the total volume of the solution to be titrated was 
made 50 ml by the addition of required volume of distilled water. 

The values of /I4, n and pL were claculated by using the Irving- Resold equations 3 ' 5 
Proton-ligand stability constant, for cephalexin was obtained by using (i) the half integral 
(from the plot of TiA vs pH) and (ii) the pointwise calculation methods; while the values 
of log K^ 1L and for the metal complexes were recorded from the formation 
curves. The refined values obtained from the least square method are given in Table 1. 

Thermodynamic parameters computed for the complexation reactions studied in this 
investigation are free energy (AG), enthalpy (Aff) and entropy (AS) changes. These have 
been calculated using the Vant Hoff isotherm, Vant Hoff equation and the Gibbs Helmhotz 
equatons respectively at 35°C and at ionic strength of O.IM with NaCICL. 

Results and Discussion 

Conductometric studies of the metal drug and metal alanine equilibria in solution 
utilising Nair and Pande’s mono-variation method 6 indicated the formation of two 
complexes in this system, with the molar ratio (metal : drug or metal : alanine) 1 : 1 and 1:2 
respectively. 


M + L = ML 
ML + L = ML 2 

where, M = the divalent metal ion, and L = cephalexin or alanine molecule or its ion. 

Charges on the species have been omitted for simplicity. Addition of third molecule 
of the ligand did not take place probably due to steric hindrance of the bulky molecules in 
accordance with the earlier reported results with penicillin and cephalosporin groups of 
drugs 7 ' 12 . 

The pH titraton curves in the eases of cephalexin and alanine lie before the held curve 
(Fig. 1), indicating that in acidic solutions each molecule of the ligand is in association 
with one equivalent of proton. This proton should be attached with the lone pair of electrons 
present on the nitrogen of the (3- lactam thiozoiidine ring in case of cephalexin and on the 
nitrogen of amino group in case of alanine as indicated in the earlier studies*'*' this 
corresponds with the two successive steps of deprotonation of the protonated ligands: 
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HLH + = H + L~ + H + 

(1) 

H + L~ = L" + H + 

(2) 


Eqn. (1) corresponds with the ionisation of the carboxyl group of the ligands, while 
equation (2) denotes deprotonation of the proton attached to the nitrogen atom. The former 
takes place in the pH range 3.5 to 5 while the later is observed between pH 6.5 to 10. 

Earlier reports 7 ' 12 also indicate similar steps of deprotonation for the protonated species 
of penicillin and cephalosporin drugs and alanine. 

Representative set of experimental titration curves obtained according to the sequence 
for different M(II) : drug : alanine systems studied reveals that below pH 6.4 formation of 
different M(II) : drug (cephalexin) binary complexes takes place. Thus the drug species, 
L* may coordinate with M(II) ions as a (N,0) donor ligand using nitrogen of OCN group 
and oxygen of carboxyl group. This is clear from the appeared divergence of each of the 
1:1 binary M(II)-drug titration curve from that of the corresponding free drug curve. 
Generally such complexes are stable upto pH 5-8 in case of different metal ions. Beyond 
these pH values, hydrolysis of the complexes leads to the formation of hydroxocomplex 
species. 

Careful examination of titration curves of the ternary complex reveals that the 1:1:1 
mixed ligand complexes are formed at pH values higher than those corresponding to the 
formation of the binary complexes of the drug (Cephalxin). This behaviour indicates that 
the binary complex ML is first formed and it coordinates with alanine(A) as secondary 
ligand to yield the mixed ligand complex MLA in a stepwise manner, similar to the formation 
of 1:2 binary complex. Further the mixed ligand curve(F) and ML(binary) curves are 
overlapped at relatively low pH values indicating that the amino acid does not bind with 
the metal ion in this range. At high pH the divergence of the two curves takes place. As 
the proton attached to the nitrogen is not removed below pH 5, hence the stepwise 
complexation equilibria can be represented as follows : 


M++ + H + L" = (ML) + H + 

( 3 ) 

(ML) + + HAH + = (MLA) + 2H + 

( 4 ) 


where H + L~ = protonated cephalexin ion, and HAH + = protonated alanine molecule. 

( a) Proton Ligand Stability Constants : The proton ligand stability constants of the 
drugs cephalexin have been calculated at 25°C and 35°C and 0.1M (NaC10 4 ) utilising the 
Irving Rossotti pH titration technique. 3,4 



metal drug equilibria in solution 
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With the increase in temperature, the values of log decreases showing an exothermic 
change i.e., with increasing temperature ligand show dissociation of the species. 

(b) Metal Ligand Formation Constants : The stability constants of 1:1 and 1:2 binary 
and 1:1:1 ternary complexes were calculated at tw o temperatures 25° and 35 ’€ and at 0.1 M 
ionic concentration, utilising Irving- Rossotti method and its modification for ternary 
systems. 3 ' 5 

The values obtained for the binary and ternary complexes are reported in Table 1. 
The ternary complexes are characterised by values of log log and Alog K vTable 

1 ). 

The results reveal that in general with the increase in the temperature log K^ iL and log 
k ml values decrease, indicating exothermic nature of complexation reactions, 

V MLA w 


M + L = ML 

(3) 

ML + A = MLA 

(4) 


A comparision of the stability constants {Table 1) of the binary and ternary complexes 
of the metal ions with the drug (cephalexin) and alanine gives interesting results. It may 
be noted that the ternary complexes, MLA, are comparatively more stable than the 
corresponding 1:2 binary complexes of the drugs. However, the ternary complexes in both 
the systems are characterised by low stability compared to the corresponding M(II) - drug 
binary complexes, i.e. K§ L > K% L U . Due to statistical, steric and electrostatic considerations 
the stability constants for the stepwise formation of 1 : 1 : 1 ternary (MLA) complex is expected 
to be smaller than the corresponding binary system and the values of Alog K 
|= log K^ l u - log should be negative 13 as obtained in the present investigation. 

According to Sigel a positive/or a value that is less negative than - 0.6 for the stability 
qualifying factor, AlogA, suggests favoured formation of ternary complex 14 . The AlogA' 
values obtained in the present study are in general less negative indicating greater 
stabilization of ternary complexes as compared to the 1:2 binary (ML 2 ) complexes of the 
drugs (Table 1). 

The order with reference to metal ions is, Mn(II) < Mg(II) < Cd(II). 
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( c) Thermodynamic Parameters : The value of changes in free energy (AG), enthalpy 
(AH) and entropy (AS) accompanying the formation of the binary and ternary complexes 
using the standard equations 15 are reported in Table 1. 

The negative values of AG show that the driving tendency of the complexation reaction 
is from left to right and the reaction tends to proceed spontaneously. 
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Abstract 

This paper is concerned with the determination of thermal stresses in different layers of a three-ia% ered 
system when the temperature is applied on the surface z = 0. This is axisymmetric, elasto-static 
problem. Transform method has been applied for the solution and numerical values are obtained for 
resultant stress. 

(Keywords : layered problem/axisymmetric/elasto-static/Hane-transform/Tieat flux/linear 
temperature distribution). 


Introduction 

This is an axisymmetric, statical three-layered problem of thermo-elasticity. This paper 
investigates thermal stresses distributed over three-layered media due to flux of heat applied 
on a surface z = 0 of the upper most layer. Layered problems have been solved by Burmister 1 , 
Nowacki 2 and double layered problem by Paria 3 due to application of concentrated force 
on the upper face of the layer. Transform technique has been taken for solution. Heat flux 
is taken to be linear. Experimental result of elastic constants are taken from international 
critical tabled 


Statement of the Problem 

Materials under consideration occupy the lower half of the plane z = 0 and the axis of 
z being taken positive w r hen drawn into the material. Physical quantities involved here all 
are symmetrical about the z axis. Lower boundaries of the first and second layeres are 
given by the plane z = h { and z = h 2 and the underlying mass extended to infinity. 
Interfaces are supposed to be perfectly rough so that stresses and displacement are 
continuous across it. 
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Method of Solutions of the Thermo-Elastic Problem 

Physical quantities involved in the problem are all symmetric about z-axis, so four 
non-zero components c r , c 0 ,a 2 ,i zr and two non-zero displacement components u, u z are 
retained in terms of stress function O satisfying the differential equation 


V 2 V 2 <J> 


\l_ Id. 31? 

dr 2 + r dr + dz 2 y 


0 = 0 


( 1 ) 


Here, stress function O due to H&nkel transform will take the form G, defined as 


Gfez) = 1^ r<I>J 0 (Z,r)dr, (E, > 0, Z > 0) 


(2) 


whose inverse is given by 

O = f 5G(5, Z)J 0 (E, r) dZ> 

0 


(3) 


Let the solution of (1) be given by 

G (q , z) = (A + Bz) (2 sin h qz + e~ + (C + Dz) (2cos h E, z - e^) (4) 


where A, B, C, D are functons of E, to be determined from the suitable boundary conditions. 

To determine the streses, the potential of thermo-elastic displacement \|/, related to the 

equations, 


3y 

3r 

3ty 

dz 


~ ( u r)r 

= ( u zh 


( 5 ) 


is considered in the steady state of temperature field given by 



THERMO-ELASTIC STRESS DISTRIBUTION 23 

V 2 r=0 (6) 

Since different layers have different thermal properties, different potential displacement 

functions are chosen. From the stress- strain relations of thermo-elasticity and the equation 
of equilibrium 3 , we have 

V 2 Vi = P,r, (7) 

1 + a, 

(3 = a,.,i = 1,2,3 

where a i are Poisson’s ratios and a i the coefficient of linear thermal expansions 
Applying Hankel transform on (6) and (7) we obtain. 

( ,2 \ 

M&z) = 0 

dz J 

and 

4 - S 2 W»z) = 

dz J 

where 

M&z) = f rTJ 0 &r) dr (10) 

0 

and 

L ; (4,z) = rv ; J 0 (qr)^r (11) 

Let the solution of (1 1) be 


(8) 


(9) 




( 12 ) 
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Then the solution of (10) is 

M&z)=A 0 e~' v . ( 13 ) 

vv here A L) is a function of q. 

Applying (12) we have, by using the eqn. (20-22) 3 

J o K<T)r Jq ( %r)dr = - + z%) A 0 e~& (14) 

J 0 r ( x zrh J i &) dr = - M/Pi (4) A o e ~^ (15) 

f Q nu r )r J x &r)dr = ± p, (1 + z^A 0 e~^ (16) 

L r ( J Or J o (&■) dr = | M A o e~ fe ( 17 ) 

Taking solution of (1) for different layers as : 

For the upper layer 

<5. :) = (A i + BiS) (2sin/t + e" 5 *) + (C, + D^) (2cos/i - e fe) (18) 
For the middle layer 

e = ( A i + B 2") (2sin/i + + (C 2 + D 2 z) (2cosA ^ - <&) (19 ) 

For //ie underlying mass 


^4. 4 = (Aj + B 3 z)e ^ £>0 , Z> hz 


( 20 ) 
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It is to be noted that stresses and displacements in the underlying mass vanish as z — > °°- 
So the components of stress and displacement for the upper and middle layers are obtained 
by using the eqn. (63-66) 5 when j=l, 2 are 

£ r (a.), J 0 (tr) dr = (l- 2ay) £ B } (2sin h & + e -5 *) + (1 - 2 a,.) q 2 D ] (2cos h & - x^) 

+ £ 3 (Aj + Bj z) (e~ ^ - 2cos/i £z) + q 3 (Cy + D } z) - 2sin/z £z) (21 ) 

£ r (x zr )j 7, (£r) dr = 2a,- £ 2 By (2cos h &-e~^) + 2a, £ 2 Dj (2sin iA §z - e^) 

+ £ 3 (Ay + Bj z) (<f & + 2sin A Cz) + S 3 ( Cy + D ; z) (2cosB qz - e^) (22) 

£ r (u z )j J x (£r) 7r = § [s,- (e“ & + 2sin/z £, z) + D } (2cos gz - <£) 

+ £ (Ay + By z) (2cosA £z - e 5 *) + l (Cy + Dj z) (2sin/i £,z - e fe )] (23) 
£ r (u,)/ 0 £r) dr = § [2(1 - 2a,) By (2 cos/icz -e*) + 2(l- 20 y) Dj (2sin/i£z - e fe ) 

- § (Ay + Byz) (2sinA £z + e~ fe )-^(Cy + Dj z) (2 cos/i^z - e fe )] (24) 

For the underlying mass 

t» ~ B 

r(a z ) 3 J Q (fy)dr = £ 3 A 3 + B 3 z + -^ (1 -2ct 3 ) e ^ (25) 

J 0 r (V)3 J \ (^ r ) dr = £ A 3 + B 3 z - ~jr B 3 e ^ 


(26) 
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J Q r ( u rh J 1 &) dr = 
J Q r ( u zh J 0&) dr = 


1+ G 3 

£3 


£3 

^ - A 3 - £ 3 z 


.-ft 


1 +c 3 


2B 

A 3 + £3 z + ~ ^ (1 — 2 ct 3 ) 



(27) 

(28) 


Boundary condition : 

In order to nullify the stresses on the boundaries the following conditions are to be 
satisfied: 


At z = 0 

- («r)r, = ( U r) 1 

“ ( U ;)t i = ( U Z ) 1 

- (a z ) Ti = (a,), 

"(V)r, = (V)i 

At Z = /ij 

“ (« r )r 2 = («r)2 

- (« z )r, = (« z ) 2 

- ( G ;) r 2 = (°;) 2 
- ( T y)r, = ( x Z rh 

At z- h 2 

~(° z ) t 3 = ( o z h J 

- (^zr)r 3 — (^jr) 3 | 


(29) 


(30) 


(31) 


Conditions (29) with the help of (14) - (17) and (21) - (24) we have, 


Aj = A,'A 0 

£1 = B/Aq 

Ci = Cj ' A 0 

Dy = D / Aq 


( 32 ) 
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where 


, 1 Pi £ i , ? s p 

A l = -Cj - 2(1 + Oi) 


B/ = D{ = 1 lit - 


2(1 + 0!) 4(1 -Ca- 


using (14), (15), (25) and (26) on the conditions of (31) we have, 
a 3 = A 3 Aq 1 


where 


b 3 = 

= B/ Bq j 


2 °3 „ 

A 3 

■f" 


H 3 P 3 

b 3 ' 



From the conditions (30) with the help of (14) - (17) and (21) - (24) we have 


A2 — A 2 Aq 
B 2 = B 2 ' Aq 

c 2 - C 2 A 0 
D 2 = D 2 Aq 


where 


A{ = ^ [j3 (1 - 2a,) - 2^A 1 } p B 2 - (1 - 2a 2 ) D 2 S 3 ] 


2(1 -a 2 )S s -g 2 5 1 

fl 2 = 2 (1 - a 2 ) (<?! -P9 2 ) 
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„ , S 2 -pB 2 (1 4o 2 + 25/jj) D 2 
- = 2' 

n ' _ Hi Z^llE^lZh H 
2 ~ 2 (1 -o 2 ) (q t -pq 2 ) 

_ d+2 cr 2 ) P 2 
1 2(1 + 02)5 

(1 + 2c/jj) jj. 2 P 2 


MjlHij _ Ei h 1 h 
\ 2(1 + 02) 5 


- \lh ~ 


EjP] M 1 + ^i) 


2(1+ o 2 )J £ 


S 5 = 2S 4 -(p- l) 2 S, 


p = e^ 


<h - (1 - 4 o 2 + 4 p- p 2 ) p 

q 2 = 2(p + 1 - 2 o 2 ) - (1 + 4a 2 ) (p - l) 2 

The relations (32), (33) and (34) give ten constants A„ £„ C„ D„ A„ C-, D A 
and B 3 m terms ofA 0 . So the formal solution of the problem is completed. ’ ’ 2 ’ 3 

The resultant stresses are in the direction of j-axis is 


On c = 0, 


= JT[(a,), + (o,) r J J 0 (fy)dr 



thermo-elastic stress distribution 


29 


At Z = /ij , 

£ r(a z )R 2 Jo&)dr = rj^ + (a z ) T ^J 0 (%r) dr 

= [- Z? A 2 e^ + Z? (I - 2a 2 - Zji x ) e^< B 2 + ? C 2 e~ ^ 

+ Z, 2 (1- 2a 2 + D 2 e' w ' - M 2 P 2 (1 + 5*i) g " y ’ 1 ] A o 


(36) 


At Z = h 7 , 


(37) 


r (c,) fij J 0 (£r) = {^(a.), + (c^r^o (^ r ) dr 

= [(1 - 2ct 3 ) £ 2 fi 3 + (A 3 + B 3 /i 2 ) Z? - Ms Ps C 1 + ^ 2 )] A ( 

Thus, the total thermo-elastic stress (a z ) R in the underlying mass is 

£ r(a z ) R J 0 &)dr = f rpa^ + (<*-)«, + (o z ) R ^ J 0 &) dr (38) 

Flux of heat on the boundaries : 

Let the flux of heat in a region of the surface a = 0, distributed through layers in the 
underlying mass be 


ar _ Jr 

dz =J \a 


0 <r<a 


= 0, r > a 

Using dimensionless variables 

£A 0 (I) = a X (&*) ’ r\ = Zfi, P = ~ 


where a is some length and tj, a new 
on z = 0 


variable of integration, we get from (10) and (13), 
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^ - f r ix 01) Jo (P^) < 39 ) 

3; a o 

the Hartkei inversion theorem 

7 (n) = - - J 1 P/.(P) J o (P 1 !) ^ (40) 

For a simple physical situation we consider a linear temperature distribution f{p) = j£p, K 

= constant, then from (40) 

With this value of %(rj), the problem is completely solved since only unknown A 0 (%) is 
now known. 

Sumerica! Results : 

If the upper layer be concrete pavement, the middle layer be gravel base and the 
underlying mass be natural soil then elastic constants for those materials are 4 

E = 2. IS x 10 s g/cm 2 , £, = 1 . 1 x 10 8 g/cm 2 , £ 3 = 0.4 x 10 s g/cm 2 

u = 5 x 10"VC, a 2 = 7.5 x 10~Vc, a 3 = 2.3 x loVc, 

u. = 0.94 x 10 s g/cm 2 , p. 2 = 0.43 x 10 s g/cm 2 , fi 3 = 0.14 X 10 8 g/cm 2 

c. =0.15, a, = 0.25 . ct 3 = 0.5 

A'. = 6.4 x 10"\ K z = 6.7 x 10' 3 , K 3 = 2.9 x 10~ 3 

So. evaluating constants for a given value of T| when a= 1, and h x =l, h 2 =A, we have 

.4. =-1471.9, £, = - 1.557, C, = - 1471.59, D, = - 1.557 

4, = - 46.43, £, = - 34.90, C 2 = - 6967.67, D 2 = - 2335.59 

.4, = - 85.56, £ 3 = - 85.56 

So applying dimensionless variables on (38) substituting the value of %(t)) from (41) and 
the values of constants in (38), we get 
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r P(a-) /o(PR)^P = - 2709.43 \ 7, 


(n) 


whose Hankel inverse transform is 


(a t ) R = - 2709.43 x K f 


J i (0) J 0 (PR) 


dr\ 


= - 2709.43 x K x F 


1 1 , vl 

2’ 2’ 1;P J' 


= - 2709.43 x Kx-, 

T C 


= - 2709.43 x K j E 


1 I. 2 .± 

2 ’ 2 ’ A p 2 

V K J 


F denotes hypergeometric function. 
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Abstract 

The thermal instability of a layer of Walters eiastico-\iscous fluid (Model ET) acted on b> a uniform 
rotation is considered. For stationary convection, Walters elastico-viscous fluid behaves like a 
Newtonian fluid. It is found that rotation has a stabilizing effect and presence of each-rotation and 
viscoelasticity introduces oscillatory modes which were non-existent in their absence. The sufficient 
condition for the non-existence of overstability is also obtained. 

(Keywords : thermal instability/ Walters elastico-viscous fluid tModei B’)/ uniform rotation,;. 

Introduction 

The thermal instability of a fluid layer with maintained temperature gradient by heating 
the underside, plays an important role in geophysics, interior of the earth, oceanography 
and the atmospheric physics etc. and has been investigated by several authors. A detailed 
account of the thermal convection in Newtonian fluids in the presence of uniform rotation 
has been discussed exhaustively by Chandrasekhar. 1 Vest and Arpaci 2 have studied the 
stability of a horizontal layer of Maxwell’s viscoelastic fluid heated from below. Bhatia 
and Steiner 3 have considered the effect of a uniform rotation on the thermal instabilty of 
a viscoelastic (Maxwell) fluid and have found that rotation has a destabilizing influence 
in contrast to the stabilizing effect on Newtonian fluid. Sharma 4 has studied the thermal 
instability of a layer of viscoelastic (Oldroydian) fluid acted on by a uniform rotation and 
found that rotation has destabilizing as well as stabilizing effects under certain conditions 
in contrast to that of a Maxwell fluid where it has a destabilizing effect. 

There are many elastico-viscous fluid that cannot be characterized by Maxwell’s 
constitutive relations or Oldroyd’s 5 constitutive relations. One such class of elastico-viscous 
fluid is Walters fluid (Model B’). Chakraborty and Sengupta 6 have studied the flow of 
unsteady viscoelastic (Walters liquid B’) conducting fluid through two porous concentric 
non-conducting infinite circular cylinders rotating with different angular velocities in the 
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presence of a uniform axial magnetic field. In another study, Sharma and Kumar have 
studied the steady flow and heat transfers of Walters fluid (Model B’) through a porous 
pipe of uniform circular cross-section with small suction. It is this class of elastico-viscous 
fluids we are interested in particularly the thermal convection in such fluids in the presence 
of rotation. 

The present paper, therefore, attempts to study the thermal instability in Walters 
elastico-viscous fluid (Model B’) in presence of uniform rotation. 

Formulation of the Problem and Dispersion Relation 

Consider an infinite horizontal layer of Walters elastico-viscous fluid (Model B’) of 
depth d which is acted on by a uniform rotation if (0, 0, Q.) and gravity force gfo, 0, -g). 
This layer is heated from below so that a steady adverse temperature gradient 
P(= \dT I dz |) is maintained. The equations of motion, continuity and heat conduction 
governing the flow of Walters elastico-viscous fluid in the presence of rotation are : 




Po 


Vp + 


v-v 


,1 

dt 


V 2 q>+ gfl + ^ 
Po 


2 (q*x if) 


V . 0 , 


§ : + (^v)r= x v 2 r, 


0 ) 

( 2 ) 

0 ) 


where q*(u, v, vv), p, p, T, v, v' and x are respectively the velocity, pressure, density, 
temperature, kinematic coefficient of viscosity, kinematic viscoelasticity and thermal 
diffusivity. The equation of state for the fluid is 

p = p 0 [l - a(T - T 0 )J (4) 

where p 0 , T 0 are respectively, the density and temperature of the fluid at the reference level 
z = 0 and (X is the coefficient of thermal expansion. 

Here we make the Boussinesq approximation under which the density changes may 
be disregarded in all terms in the equations of motion except the one in the external force. 



THERMAL INSTABILITY IN WALTERS ELASTICO-VISCOUS FLUID 


35 


The initial state is one in which the velocity, density, pressure and temperature at any pont 
in the fluid are respectively given by 


Y= (0,0,0), p = p(z),p=p(z)tand T = T(z) 


(5) 


Let q*(u, v, w) 8p, 8p and 0 denote respectively the perturbation in velocity (0,0,0), 
density p, pressure p and temperature T. Then the linearized perturbation equations are 


da* 1 „ . 

-? = - — V dp + 

dt p 0 


V - V' 


, l 

dt 


V 2 q*+ + 2(?x Qf, 

Po 


( 6 ) 


v . r= o , 


(7) 


| + (?V)T=xV 2 0. 


( 8 ) 


Within the framework of Boussinesq approximation, eqn(6)-(8) becomes 

d v 2 w 


dt 


- gcc 


fd 2 e a 2 e) __ac f , d' 

[dx 2 dy 2 ) dz [ dt 


V 4 w 


(9) 


f -2£2^ 

dt dz 


v - v 


, A 

dt 


v 2 ; 


( 10 ) 


— - y V 2 
dt X 


0 = P W , 


(ID 


where V 2 = 


"a 2 a 2 a 2 ^ 

dx 2 + dy 2 + dz 2 


and C ^ f~ “ denotes the z-component of vorticity. 


We now analyze the disturbances into normal modes, assuming that the perturbation 
quantities are of the form 


[w, 0, Q = [W(z), 0 (z), Z(z)l exp (ik x x + ik y y + nt ) 


( 12 ) 
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where k x , A; are wave numbers along w and y directions respectively, A 2 = + A 2 and n 

is, in general, a complex constant. 

Letting 


, , nd 2 v 

a = kd, a = — , p , = - , 


(x, y, 


z) ( xd, yd , zO 



Equations (9)-( 11) using expression (12) give 
ct (D 2 - a 2 ) IV + 




a v y 


rl/2 
•A ' 


[D- - a 2 - p . a] 0 = - 


z = 

(1 - Fa)(D 2 - a 2 ) W , 

(13) 


If 



f DW , 

a 

(14) 

p/] 

X 

w , 

(15) 


where T A = 


4a 2 / 


•> denotes the Taylor number and D = — , 
v z dz 


Applying the operator [D~ — a" - py o] [(1 — Fa) (D 2 — a 2 ) — cr] to eqn. (13), we 
eliminate © and z to obtain 


a (D 2 - a 2 ) (D 2 - a 2 - Pl a) [(1 - Fa) (D 2 - a 2 ) - a] W 
- R a 2 [(1 - Fa)(D 2 - a 2 ) - a] IV - T A (D 2 - a 2 - Pl a) D 2 W = (1 - Fa) 
(Ft - a 2 - p, a) (D 2 - a 2 ) 2 |](l - Fa) (D 2 - a 2 ) - cr] IV , (16) 
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where R 


g a M 

VX 


stands for Rayleigh number. 


We now assume that the fluid layer is confined between two free boundaries. The case 
is of artificial nature but due to mathematical simplicity it enables us to show the behaviour 
of rotation on Walters elastico-viscous fluid analytically. Then the boundary' conditons 
appropriate for the problem are 


W = d 2 W = 0’, 0 = -0\ DZ = 0 at z = <0 and z = 1 . (17) 


Dropping the caps for convenience and using the above boundary conditions, it can be 
shown that all the even order derivatives of W must vanish on the boundaries and hence 
the proper solution of eqn. (16) characterizing the lowest mode is 


W = W 0 sin 7i z , (18) 

where W 0 is a constant. 

Substituting (18) in eqn. (16), we obtain the dispersion rotation 

(1 +x)(i +x + i,p l a,) ((1 - 4 crj jc 2 /0 (1 +*) + c<L) 

((1 -i <7, TC 2 F) (1 +*) + 4C 1 ) + 7’,(l +x + i Pl a,) 

R _ (19) 

JC [(1 - i CTi It 2 O (!+■*) + 6 ®il 


where 


«i 





it being remembered that a can be complex. 

Here we consider the overstable mode and so c l is real in eqn. (19). 

For the case of stationary convection i.e. 0 = 0. eqn. (19) reduces to 
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(1 + *) 3 + T x 


( 20 ) 


[a result given by Chandrashekhar 1 . eqn(130), p. 95]. 

We thus find that for the stationary convection, the visco- elasticity parameter F vanishes 
with a and Waiters elastico-viscous fluid behaves like an ordinary Newtonian fluid. 

To study the effect of rotation, we examine the nature of dR l /dT v From eqn. (19), it 

follows that 

_ (1 + x + i pi a,) 

dT \ ~ *[(1 - +*) + iaj (21) 

On equating real and imaginary parts of eqn. (21), it follows that 
dR x i 

dT\ = x ’ ( 22 ) 

which is always positive. The rotation thus has a stabilizing influence. 

Stability of the System and Oscillatory Modes 

Multiplying eqn. (13) by W*, the complex conjugate of W, integrating over the range 

of z and making use of eqn. (14) and (15) together with the boundary conditions (17) we 

obtain 


- al. + 


sm: 

vp 


(/ 2 + p l a* / 3 ) - d\ 1 - Fo*) / 4 - d 2 a 1 5 = (1 - Fa) I a , 


(23) 


h = (|DSf + o ! |S| 2 )* , 


where 
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/ 3 = Jj©i 2 *. i4 = f Q {m 2 + a 2 \z\ 2 y z , 

h = 1* M* dz , / 6 = J* (|D 2 W| 2 + 2 a 2 |DWj 2 + a 4 |W| 2 ) dz , (24) 


and a* is the complex conjugate of a. The integrals l x - I a are all positive definite. Putting 
a = a r + i cr, in eqn. (23) and equating the real and imaginary parts, we obtain 





+ 


vP 


Pl / 3 + d 2 FI 4 - d% + FI 6 





(25) 


and 




h + 


me * 

vP 


Pl / 3 + d 2 FI 4 - d 2 / 5 - F/ 6 


= 0 , 


(26) 


It is evident from eqn. (25) that o r is positive or negative. The system is therefore stable 
or unstable. It is clear from eqn. (26) that rotation and viscoelasticity introduce oscillatory 
modes in the system which were non-existent in their absence. 


The Case of Overstability 


Here we discuss the possibility of whether instability may occur as an overstability. 
When the marginal state is oscillatory, we must have o r = 0. a ; * 0. 

Since for overstability, we wish to determine the critical Rayleigh number for the onset 
of instability via a state of pure oscillations, it will suffice to find conditions for which eqn. 
(19) will admit of solutions with a r real. 

Separating the real and imaginary parts of eqn. (19), we have 


R x (a - 1) = <? - aaj [(n 2 Fa - l)(n 2 Fa ~ 1 - 2 + T l , (27) 

R x (a - 1) (1 - n 2 F a) = (1 - n 2 Fa) a [2a 2 - a?p,(l - n 2 Fa)] 


+ Pi (a 3 + T { ) , 


(28) 
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where 1 + .v = cc. 

Eliminating Rj between eqn. (27) and (28), we have 

[(1 - r 1 2 Fa) 2 a (1 + P[ -n 2 Fa)]a\ + [a 3 (1 + P{ - n 2 F a) - T x (1 - Pl - n 2 Fa) = 0. (29) 
Eqn. (29) is of the form 

A a] + B = 0 . (30) 


where A,B are coefficient of a] and constant term respectively in eqn. (29). As a, is real 
for overstability, the values of q is positive. Eqn. (30) shows that this is impossible if A 
> 0 and B > 0. Therefore, A > 0 and B > 0 gives sufficient conditions for the non-existence 
of overstability, which yields 


Pl > it 2 Fa - 1 


and Pl < 1 - n 2 F a 


(31) 


which further implies that 


F(k 2 + k 2 d 2 ) - 


< - < 1 - F (n 2 + k 2 d 2 ) 

1 


(32) 


Thus for F(tc 2 + k 2 d 2 ) - 1 



1 - F (n 2 + k 2 d 2 ), overstability cannot 


occur. The condition (32) is, therefore, sufficient condition for the non-existence of 
overstability, the violation of which does not necessarily imply occurrence of overstability. 
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Abstract 

Within the framework of two parameter octahedral potential for a dipolar impurity in ionic crystals, 
it is possible for impurity to possess minimum-energy orientational configurations in two 
crystallographic directions simultaneously. Two possible combinations are <1 1 1>— <1 10>and<l 1 1> 

- <100>. Taking second model here we have tried to findout the expressins for polarization. 

(Keywords : ionic crystals/octahedral potential/polarization) 

Introduction 

Diatomic and off-centred monoatomic impurities in crystals with octahedral symmetry' 
have attracted considerable attention. Within limit of infinite potential barriers these 
impurities align themselves along certain preferred crystallographic directions. In case of 
finite potential barriers, they tunnel between the equilibrium orientional configurations 
resulting in interesting electrical, thermal and optical properties. The above mentioned 
properties, depend critically on the actual minimum - energy orientational configurations 
of the impurity, apart from certain other parameters such as the height of the potential 
barriers etc. The simplest form of the octahedral potential for the angular motion of diatomic 
impurity in solid-state matrices was first presented by Devonshire 1 and is given by : 

V(0, a) = - (5m) [(/ + / + z 4 ) / r 4 ] (1) 

where xlr — sin0 cosoc, y/r = sin0 sinoc and z/r = cos9. Many experimental results for 
systems such as OH' and CN* ions in alkali halide matrices 2,3 and impurities like HC1, 
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HBr etc. in rare gas matrices 4,5 have successfully intepreted in terms of this one-parameter 
Devonshire potential. It is easy to show that for this potential the minimum-energy 
orientational configuration for the impurity are along one of six <100> directions and 
along one of the eight <1 1 1> directions if K is positive and negative respectively. A <1 10> 
minimum-energy configuration was not possible within the framework of the Devonshire 
potential This potential in fact, retains only the first significant term (the L = 4 term) in 
the expansion of octahedral potential. To make it for the <110> configuration, the next 
term (the L=6 term) was added to make the two parameter octahedral potential as given 
by 


m a) = - (5K/2) [(/ + / + z 4 )/r 4 ] + (l/2)/d [180 (x 1 y 2 z 2 /r b ) - (15/r 6 ) 

(x 2 y 4 + x 2 z 4 + y 2 z 4 + z 2 x 4 + z 2 y 4 ) + (2/r 6 ) (x 6 + / + z 6 )] (2) 

It was concluded in a previous communcation 6 that if (10/21) < (jd/K) < (-20/189), a 
<110> equilibrium orientational configuration is physically possible. This was also 
subscribed later by Beyeler 7 whose condition was 26° <a< 173°, where a = tan' 1 (tf/K). 

The nature of Li + in KC1 matrix are explained by single tunneling model proposed by 
Gomez et al. 9 but not well applied for other’s as NaBr:F~ and RbCl:CN“ systems. Various 
attempts 0, to explain the experimental data were made using <1 10> model, but got limited 
success. Experimental results of RbChCN system were not explained using this model. 
This system was also not very well explained in terms of the <110> and <11 1> model 
proposed by Pandey et al. 1 In the present paper we have attempted to find out the expression 
of polarization of < 1 00> and < 1 1 1 > model. 

Limitations of Old Models and Scope of Simultaneous 
Potential Minima in two Crystallographic Directions 


Inspite of KCL:Li + system, which shows agreement with the <1 1 1> off-centred model 
of Gomez et al. , many other paraelectric systems present difficulty in understanding the 
various experiments associated with them. In those cases existence of two sets of 
minimum energy orientational configurations of the impurity in the matrix give interesting 
and exciting results. Initial experiments of Seward and Narayanmurti 3 presented evidence 
that lmpunty could have potential minima along any of the six <100> directions. Using 
the model RbChCN system, paraelectric resonance 10 and specific heat 11 at low temperature 
could not be explained. However, Luty 12 and Beyeler 13 on stress-induced splitting of the 
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Fig. 1 - Pocket states for simultaneous potential minima along the <1 1 1> and <100> directions in the Oh 
symmetry group. 

intramolecular vibrational line showed that correct minimum-energy orienta- 
tional configuration of CN" impurity in ionic crystals was along one of eight <11 1> 
directions. The model fails to explain the double bump in specific heat results of RbCkCN' 
system. 

Similarly, the calculation of Quigley and Das 14 shows minima position along <11 1> 
direction while the high field polarization and paraelectric cooling experiments of the 
Rollefson 15 are better fitted along the <110> direction. After some time Roliefson 16 
suggested the existence of potential minima along the <11 1> direction, which are only 
slightly less deep than the <1 1G> minima. This thing is physically possible even within the 
framework of the two-parameter potential of eqn. (2) discussed above. Specific heat 
anomaly at low temperatures was not fully explained within the simple <1 10> model 
There are several examples which could not be fully explained by simple tunneling model. 
At very low temperatures, when most of the molecules remain in lower well, the 
behaviour of system can easily be explained in terms of simple model. But at elevated 
temperatures, more and more molecules are transferred to the upper well and the system 
shows a marked departure. In some cases when both wells are equally deep, the behaviour 
of the system becomes crucially linked to the experimental conditions and one gets 
controversial results. 
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Theory 

There were six equivalent <100>, the eight <1 1 1> and twelve <1 10> states in paper 
of Gomez et al. 9 In present case there is a set of each of six and eight equivalent states for 
the case of the simultaneous occurrence of the potential minima along the <100> and <1 1 1> 
directions. The designation of states is given in Fig. 1. 


For a finite crystal field, there will be matrix elements of the Hamiltonian connecting 
the different states i.e., there will be a finite probability of tunneling through the potential 
barrier to a new orientation. There will be matrix element connecting the states, for example, 
\A> of the <1 1 1> minima group to the nearest neighbouring wells |B>, |D> and |£> of the 
same minima group. Similarly, there will be equal matrix element connecting the states 
|a> of the <100> minima group to the states \c>, \d>, \e> and |/> of the same group. Here 
we assume that the tunneling matrix elements of the far off well are insignificant compared 
to the tunneling matrix elements to the nearest-neighbouring wells. Now impurity can 
tunnel from a well of the <1 1 1 > group of the nearest-neighbour well of the <100> group. 
Thus we have equal matrix element of the Hamiltonian connecting the states [A> of the 
<1 1 1> group of the states |< 2 >, |c> and \e> of the <100> group. 

Explanation of Field Polarization Effect 

If we apply an external field E at an angle 0 with the dipole orientation, then the average 
field polarization (P) in the field direction is given by 


n 


P = N\i 


£ cos9 r e ^ EcosQ r /kT 

r- 1 


^iEcosQ r /kT 


r=l 


( 3 ) 


whe« i - Boltzmann Constant, T - Absolute temperature, M = Dipole concentration, n = 

Dipole moment. 


Ut ■* = ■ So eqn ‘ (3) takes the 
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Table 1 - Valuesofcos0randcorrespondingdipolarenergyfor<lOO>systemwithappliede!ectricfieldin<lOO> 

direction. 


Well 



010 

010 



cos 0r 

1 

-1 

0 

0 

0 

0 

-p£ cos 0r 


\iE 

0 

0 

0 

0 


Table 2 - Values of cos 0r and corresponding dipolar energy for<l 00> system with applied electric field in <1 1G> 
direction. 


Wells 

100 

100 

010 

010 

001 

001 

cos 0r 

l/C 

-i/c 

l/C 

-l/C 

0 

0 

-p£ cos 0r 

-p£C 

p£/C 

-pE/C 

p£C 

0 

0 


Table 3 - Values of cos Grand corresponding dipolar energy for<100> system with applied electric field in <1 1 1> 
direction. 


Wells 

100 

100 

010 

010 

001 

001 

cos 0r 

-p£ cos 0r 

l/C 

-p£/C 

-l/C 

p£/C 

l/C 

-pE/C 

-l/C 

p£ /C 

l/C 

-p£/C 

-l/C 

p£/C 


Table 4 - Values of cos 0r and corresponding dipolar energy for <100>+<1 1 1> system with applied electric field 
in <100> direction. 


Wells 

100 

100 

010 

010 

001 

001 

111 

111 

111 

111 

111 

111 

ill 

111 

cos0r 

1 

-1 

0 

0 

0 

0 

l/C 

l/C 

l/C 

-1/C 

l/C 

-l/C 

-l/C 

-i/C 

-p£cos0r 

-HE 

pE 

0 

0 

0 

0 

-p£/ 3 

-p£/3 

-pE3 

p£/3 

— pt£/3 

ji£/3 

pE/3 

pE3 


Table 5 - Values of cos 0r and corresponding dipolar energy for <1 00>+<1 1 1> system with applied electric field 
in < 1 00> direction. 

Wells 100 100 010 01 0 001 001 111 HI HI HI 111 111 111 1H 

cos0r l/C -l/C l/C -l/C 0 0 2/C 2/C 0 0 0 0 -2/C -2/C 

-pEcosSr 2 pE/V 2 -pE/C p£/C 0 0 -2pE/C2p£/C 0 0 0 0 2p£/C 2p£/C 
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Table 6 - Values of cos Brand corresponding dipolar energy for <100>+<1 1 1> system with applied electric field 


in <11 1> direction. 


Wells 

100 100 010 010 001 001 

111 

111 

111 

111 111 

111 

111 

111 

cosBr 

iNT-i/VI l/VT -ihfi u€ -l/VT 

1 

1/3 

1/3 

1/3 -1/3 

-1/3 

-1/3 

-i 

-p£cos9r 

-y£ly$E/- ■fi-nEI'H \ lEI'IT -hE/V 3 liE/'JT 

-V£ 

-p£/3 -p£/3 -p£/3jx£/3 

|a£/3 

|i£/3 



Table 7 - Expressions for polarization for the simple <100> off-centred model and the present model of 
simultaneous potential minima along the <100> and <1 1 1> directions. Here * = \xE/kT and the asterisk 
denotes experssions which represent saturation value for the case x > 1 (i.e., for larger field). 

Direction Polarization as Polarization as expected from present 

of applied expected from model 

electric the simple <100> 

field model 


<100> 


P = N 


x -x 
e - e 


(Xl (/I - r 1 ') + V3 n 2 (e x ^~ e~ x ^) 
P ~ N 4 + /> + «■*!+ 3 </*** + f x ^) 


< 110 > 


< 1 ! 1 > 


P = Ny.* 


p=Nvi ; 


x/'l2_ -xl<2 
1 +e +e 



_ xtiJ _ -x/^3 
e -he 



w faiN IXe^-e x ^)+( 'l2/3)n 2 (e 2x ^-e 2x ^) 

P = N 3+ /i N2 + ^> /V2 + e 2r ^ 6 +e -^ 

P = ^L 
2 

(e x ^- e~ x ^) + ^ (/2 _ e~ x 2 + 

3(/ ]/ ^+ e - * 1 ^ 3 ) + e x ^+e~ X2 + 3(e x ^ 3 + 7^) 

p = n£ 


X COS0 r ^ cos0 r 

p = Mi ^ (4) 

X ^ cos6 r 

r= 1 


Now using eqn. (4) we find expression for polarization P in different cases, 
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Case 1 : System <100> and electric field along <100> direction : For <100> system 
with electric field in <100> direction value of cosG^ and their dipolar energy is given in 
Table 1. So using Table 1 and eqn. (4) we have expression for polarization given as 


P 


N 


e x - 

4 + e x + e~ x 


( 5 ) 


For x » 1 (i.e., » kT) 

P s ~ ( 6 ) 

where P s is saturation polarization. 

Case 2 : System <100> and electric field along <110> direction : Using Table 2 and 
eqn. (4) we have expression for polarization given as 


P 


N\i 


2 [1/V2) (e*® - 1 

2 + 2(e^ 2 + e'^ 2 ) 


p = ApyV 2 


e xti2 _ e -xtd 

1 + 7^ + e'*® 


Forx»l (i.e., \iE»kT) 

P s = N\ii<2 


( 7 ) 


( 8 ) 


Case 3 : System <100> and electric field along <111> direction : Using Table 3 and 
eqn. (4) we have expression for polarization given as 


P = N\i 


- e~^)} 

3(P* + e"^ 3 ) 


_ _ e ~*& 
-g 

e + e 


( 9 ) 
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For .v » 1 (i.e., p£ » kT) 


P s = NyJ'fi 


( 10 ) 


Now we have extended our calculation to find out the expressions for polarization for 
system <100> + <1 1 1> (present model). 

Let p and p, be the dipole moment for <100> and <1 1 1> wells respectively. 

Case 4 : System <100> + <111> and electric field along <100> direction : Using 
Table 4 and eqn. (4) we have expression for polarization given as 

p, («*< - «"*>) + 3 {(p 2 / VT) (e x ^ - e^H 
P ~ N 4 + e X| + e~ x > + 3(e^ 3 + e' x ^ 3 ) 


Pj (e x ' - e *>) + V3* pi 2 (e x ^ - e x ^) 
P N 4 + ^‘ + <T X > + 3 (e x ^ + e~ x ^ 3 ) 


(ID 


P s = ^ 


( 12 ) 


Case 5 : System <100> + <111> and electric field along <110> direction : Using 
Table 5 and eqn. (4) we have expression for polarization given as 


2|(p 1 /V2)(<? I < /VI - 

e" x ' /2 )J + 2 )(2p2/V6) (e^ - e _2x ^) 


2 + 2(e x > ,N ' 2 + 

e^ 2 ) + 4 + 2(e 2x ^ 6 + e'^) 


(p,/V2) (e x ^ - e 

~ x ^) + (V2/3) pj (e 2 ^ - 

(13) 

3+T^ 2 

+ + £ 2X 2 / ^ 6 + e ~ X ^ 6 

P s = NyJ2 


(14) 


Case 6 : System <100> + <11 1> and electric field along <U1> direction : Using 
Table 6 and eqn. (4) we have expression for polarization given as 
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3 (|X 1 /V3)’ (e x>/ ^ - e x ' 1 ^) + ^ {(«* -e x *) + (3/3) (e'^' 3 - e"^ 3 )) 
P ^ 3(^' 7 ^ + e _Jt ‘ / ^ 3 ) + ^ + e~ x - + 3 (e^ 3 + e~ x ^) 

V3" jij (^i 7 ^ 3 - e~ x '^) + ji, - e _jr - + W 3 - e - ^ 3 
^ ^ 3(e x ^ + e~ x '^) + e* + f x i + 3 e (/T 3 + e - */ 3 ) 


(16) 

Now expressions for polarization for <100> and <100> + <11 1> model (present model) 
are listed in Table 7. 


Results and Discussion 

The expression for polarization for the <100> off-centred model for impurity is 
summarized in Table 7. The same table also lists the expressions for polarization for present 
model of simultaneous minima along the <100> and <1 1 1> directions. The off- centre 
displacement along the <100> and <1 1 1> directions, the saturation polarization for fields 
along the <100> and <11 1> directions should be equal (cf. Teble 7 for the case jij = JT). 
The larger polarization for the <100> field in this case only indicates that the off-centre 
displacement of impurity of the <100> minima direction is larger than that in the <1 1 1> 
direction. 
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Abstract 

Molecular interactions between a pair of molecules of the homologous series para-n-alkoxybenzoic 
acid have been evaluated using quantum mechanical methods. Both, stacking as well as in-plane 
interacting configurations have been taken into account. Results obtained have been used to discuss 
the molecular forces responsible for the mesogenic behaviour of the systems. An attempt has been 
made to explain mesogenic and non-mesogenic characters exhibited by different members of the 
same homologous series in terms of the intermolecular interactions. 

(Keywords ; liquid crystals/molecular interactions/CNDO/mesogen/alkoxybenzoic acids). 


Introduction 

Molecular ordering in mesomorphic substances has been receiving attention for quite 
sometime 1 . A number of theoretical models to explain the liquid crystallinity have been 
reported from time to time 2 ’ 5 . These models, largely concerned with nematogenic materials, 
attempt to explain the observed orientational order for such systems. The molecular field 
approach 4 presently deemed to be the most popular and successful one is based on the 
following 6 : 

(i) The assumption of an induced-dipole induced-dipole interaction betw een molecules 
of the liquids as the sole cause of transition, 

(ii) The assumption that the probability of any configuration of centres of mass is not 
affected by the interaction between different molecular orientations, and 

(iii) The molecular field approximation to treat the resulting problem of interacting 
orientations. 
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Later modifications have included parameters representing repulsive or polarization 
potential which were essentially qualitative in nature 7 ' 8 . Commonly, a simple anisotropic 
pair potential, U(r l2 , 0 I2 ) = - U(r n ) P 2 ( cos 0 12 ) is used; where U is a scalar potential 
P'S cos e, 2 ) is a second Legendre polynomial, r 12 and 0 12 are the separation and angle 
between the molecules 1 and 2 respecively 9 . Also, the cylindrical symmetry of the rod like 
molecules is assumed together with the equivalence of 0 and (n - 0). These assumptions 
based on some experimental evidences, help to reduce computational complexities. 


The peculiar changes-characteristics of mesomorphic behaviour which occur at phase 
transitions, are primarily governed by the intermolecular interactions acting between sides 
planes and ends of a pair of molecules'. In view of the key role of molecular interactions 
m mesogenic compounds, semi-empirical calculations have been emphasized and carried 
out by several workers with an aim to explain liquid crystallinity 10 ' 12 . Tokita et al. n used 
Lennard-Jones potential to evaluate intermolecular interactions between a pair of pure 
nematogens and attempted to correlate their results with those of the molecular field theory 5 
However, rt has been observed that ’6-exp’ type of potential functions are found to be more 
ettective in explaining the molecular packing instead of Lennard-Jones potentials 13 It is 
therefore, considered worthwile to examine the anisotropy of interactions between a pair 
o mesomorphic molecules. In our earlier communications, using Rayleigh-Schrodinger 

fntermnl ti0n i ^ • m0dlfied With muld centered-multipole expansion method, 
intermolecular interaction energy has been evaluated in case of several mesogens 14 ' 18 Also 

~ made r° COrT ? t£ thC meSOgenic P r °P erties with molecular interactions! 
S “ d,eS a molecular i“' of Mowing systems have been reported in this paper 

(i) para-n-methoxybenzoic acid (PMB) 


fii) para-n-propoxybenzoic acid (PPB) 


(iii) para-n-octyloxybenzoic acid (POB) 


ring JS- ‘ir m0leC " leS is th “ Ihe )' P° sse ss a polarizable atomalie 

one e.d“I^S <7 T * *“*> “ 

head to head nai™ dim ^ u a * Th*** aCld molecuIes are capable of forming 

molecules 19 . An examin a tin he fih yd n 0Sen b ° ndS betWeen the “ C00H end groups of the 

that PMB is a pure liouid wh° h 5 C erm ° dynamic data for the homologous series indicates 
pure liquid which gives a clear melt at 184’C for C -> I transition 20 . PPB 


transforms to a nematic mesophase exhibiting C- 


121*C 


— > c, c- 


I transition 2 
146.5°C 


-» N, and 
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' > I transitions while POB shows both, nematic and smectic mesophases at different 

temperatues as under : 

75 °C 101 o C 108°C 146°C 

C >C, C »S, S »N, N *1 

Theoretical Method 

The choice of an adequate method for evaluation of interactions between a molecular 
pair seems to be important. The total intermolecular interaction energy (E tot ), neglecting 
the resonance and magnetic interactions, between a pair of molecules can be expressed as: 

Etc = E el + E^, + E disp + E rep ( 1 > 

where E el , E pol , E disp and E rep represent electrostatic, polarization, dispersion and repulsion 
components respectively. 

According to Rayleigh-Schrodinger perturbation treatment for long-range interactions, 
the electrostatic term is obtained as a first order perturbation energy while the induction 
and dispersion energy components are obtained as the second order perturbation terms" 
Repulsion energy is calculated separately either from the repulsive part of Lennard-Jones 
or Buckingham type of potentials or from an optimized potential function including the 
repulsive and attractive (dispersion) terms 23 . 

The aforesaid method has difficulties because the complete process of evaluation of 
energy has to be carried out separately for each pair of molecules while a more practical 
method would be to develop certain parameters in terms of individual molecular properties 
such as multipole moments, polarizability etc. which may reduce the conceptual 
complexities and computational complications. In literature, a number of methods based 
on R-S treatment are available on these lines 24 . Earlier workers ^employed one centre 
expansion of the multipole moments to evaluate electrostatic energy " . This method was 
found to be suitable for smaller molecules while for larger systems, however, intermolecular 
separation can no longer be treated very large as compared to the molecular dimensions. 
Also R-S treatment has some inherent inadequacies at short ranges as discussed by 
Claverie 27 ; 

An alternative approach, commonly referred to as multicentred- multipole expansion 
method appears to be more appropriate for larger molecular systems" . In this approach, 
instead of taking the centre of mass of the whole molecule as the centre of expansion. 
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multipoles corresponding to each atomic centre of the molecule, are evaluated using all 
valence molecular orbital methods such as SCFcr + Huckel 7t, Del Recr + Huckebt, CNDO, 
EHT, IEHT etc. However, CNDO/2 method is preferred because of its simplicity and 
capability to predict correct dipole moment etc 29 . Usually, the evaluation of atomic net 
charges and dipoles is sufficient to treat the interaction problem with a reasonable accuracy. 
A number of successful applications to biomolecular interactions have placed the method 
on a strong ground to account for the short and intermediate range interactions of large 
molecules 30 ' 36 . In the present computation, the method developed by Claverie has been 
followed. 

Electrostatic energy : According to the multicentred-multipole expansion, the 
electrostatic term is expressed as 

E eI = E qq + E qmi + E n , mi + (2) 

where E qq , E qmi and E mimi etc. represent monopole-monopole, monopole-dipole, 
dipole-dipole and other higher order interaction energy terms respectively. The electrostatic 
component arising due to higher order multipole moments also exist but only an evaluation 
of electrostatic energy upto the first three terms has been found suitable for most of the 

cases 37 . 

Again, monopole-monopole energy is given as : 


v 



(3) 


where q t , are monopoles at each of the atomic centres of the interacting molecules while 
r t) is the inter-atomic distance. The constant, C, is a conversion factor, approximately equal 
to 332, which expresses the energy in kilocalories per mole of the dimer. 

The monopole-dipole energy is expressed as : 

Eqmi ~ C Yj Qi P/ ■ ~3 W 

i.j ' 


while dipole-dipole term is given by 


U r 


PS? - 3 


p 9 .- 

V J 


U r* 

^- r 

\ M 


(5) 



LIQUID CRYSTALLINITY OF ALKOXYBENZOIC ACIDS 


57 


where p.,, Pj represent the atomic dipoles, the subscript of r has been removed without any 
change in its meaning and other notations have the same meaning as in eqn. (3). 

Polarization energy : The polarization energy of a molecule, say ( s ), is obtained as a 
sum of the polarization energies for the various bonds 

(s) 

£$ = C(-1/2)£ k e^AMeM (6) 

where b“ = £ M 7 4- (7) 

t *s X 

is the electric field created at bond u by all the surrounding molecules and A is the 

^ u 

polarizabilty tensor of this bond. is the vector joining the atom X in molecule (r) to the 

‘centre of polarizable charge’ on the bond u of molecule (j). 

Dispersion and repulsion energy : Dispersion and repulsion terms are calculated 
together using Kitaigorodskii 38-39 type of formulae 

E dlsp + Erep = I' 0 E (2> E(k.V) (8) 

X v 

here E (X, v) = K x K v (- A/z 6 + Be" 17 -) 

and z = V4 1 K = V(2^) (2J® (9) 

Here R% and are van der Waal’s radii of atoms X and v respectively. The parameters 
A, B and y do not depend on the atomic species, but R\ v and the factors K x and A' v depend 
on the atomic species involved. The recent values of these parameters have been taken 
from literature 40-41 . 

The molecular geometry of PMB has been taken from literature 4 ". PPB and POB have 
been constructed using crystallographic data together with the standard values of bond 
lengths and bond angles. Complete neglect of differential overlap method, (CNDO/2), has 
been employed to compute the net atomic charge and dipoles on the various atomic centres 
of the molecules 43 . 

The energy minimization has been earned out for both stacked and in-plane interactions 
separately in all the systems under study. One of the interacting molecules is kept fixed 
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16 is 



Fig. 1 - Molecular geometries of systems studied alongwith various atomic index numbers, (a) 
p-n-methoxybenzoic acid (PMB), (b) p-n- propoxybenzoic acid (PPB), (c) p-n-octyloxybenzoic acid 
(POB). 



F*g. - \ arsation of stacking energy with respect to interplanar separation. Energy is expressed in kcal/mole 

o imer. (a) p-n- methoxybenzoic acid, (b) p-n-propoxybenzoic acid, (c) p-n- octyloxybenzoic acid. 
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iroughout the process while variation in both, angular and lateral positions are given in 
ie other molecule in all respects relative to the fixed one. Accuracies upto 0.1 A° during 
liding (or translation) and 1° in rotation, have been achieved. All the computations were 
arried out on CDC ‘Cyber - 170’ computer at Tata Institute of Fundamental Research, 
dumbai. 


Results and Discussion 


Molecular geometries of PMB, PPB and POB alongwith their atomic index numbers 
ire shown in Fig. 1 while the net atomic charge and dipole moment components 
mrresponding to each atomic centre are listed in Tables (1-3) respectively for PMB, PPB 
ind POB molecules. The total energy, binding energy, total dipole moment and its 
components are given in Table 4 for all the homologues under study. 

As is evident from Tables (1-3), the oxygen atoms of carboxylic groups O 10 of methoxy 
and propoxybenzoic acids and 0 9 of octyloxybenzoic acid molecule have a constant value 
for the net atomic charge (- 0.35 e) while the other oxygen atom O n of PMB and PPB 
and O 10 of POB also have a constant charge equal to -0.225 e. Since the carboxylic oxygen 
atoms' are associated with greater negative charges and attached to one end of the 
p-n-aikoxybenzoic acid molecules, they are always capable of forming hydrogenbonded 
dimers whenever any proton is available to them - a fact observed experimentally '. Also, 
the total dipole moment for benzoic acid molecules (Table 4) shows a fairly constant value 
while the magnitude of the total energy has the following order : 

PMB < PPB < POB 


This indicates that addition of methylene group confers additional stabdity to the 
packing of the members of this homologous series. The interaction energy between a pair 
of molecules has been analysed in terms of several parameters as described below. 


Separation: This has been carried out by keeping the long molecular axes of the 

interacting molecules parallel to each other. In this case, the coord, nate axes (XX 3 are 

defined such dtat X-axis is directed along to long molecular axts wttlt K-ax.s being 

perpendicular to it and lying in the molecular plane while Z-tms is ^ndtcular m to 

lane of the molecule and passes through the centre of gravity of the molecule, which ,s 

chosen as the origin. To obtain the most stable configuration (maximum attracnon 
cnosen as me o g stacked or (b) m-plane interacting 

condition), the molecules are considered either m w w K + 

J £ . , • „ jntoroptions the interplanar separation between a 

configurations. In the case of stacking interactions, me uuc p f 

. , r nn wards at an interval of 0.5 A and at each point 

pair of molecules is varied from 3.5 A onwaras <u 
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Table 1 - Molecular charge distribution of p-n-methoxybenzoic acid. 


Atom No. Symbol 

Charge 

(e.u.)* 

Atomic dipole components Debye) 

x y 1 

I 

C 

0.392 

-0.217 

-0.221 

-0.014 

2 

c 

- 0.063 

0.070 

-0.034 

-0.049 

3 

c 

0.055 

0.060 

-0.146 

-0.001 

4 

c 

- 0.061 

- 0.084 

-0.170 

0.005 

5 

c 

0.210 

0.232 

0.033 

- 0.027 

6 

c 

-0.068 

-0.057 

0.157 

0.016 

7 

c 

0.047 

0.080 

0.121 

-0.001 

8 

0 

-0.210 

- 0.626 

- 1.210 

-0.069 

9 

c 

0.122 

-0.198 

0.244 

0.034 

10 

0 

-0.351 

0.644 

- 1.209 

0.029 

11 

0 

-0.255 

-0.555 

1.157 

0.087 

12 

H 

0.167 

0.000 

0.000 

0.000 

13 

H 

0.004 

0.000 

0.000 

0.000 

14 

H 

0.015 

0.000 

0.000 

0.000 

15 

H 

0.006 

0.000 

0.000 

0.000 

16 

H 

-0.001 

0.000 

0.000 

0.000 

17 

H 

-0.006 

0.000 

0.000 

0.000 

18 

H 

0.003 

0.000 

0.000 

0.000 

19 

H 

-0.006 

0.000 

0.000 

0.000 


*e.u, = electron unit 
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Table 2 - Molecular charge distribution of p-n-propoxybenzoic acid. 


Atom No. 

Symbol 

Charge 

Atomic dipole components (Debye) 



(e.u.)* 

X 

Y 

z 

1 

C 

0.393 

- 0.220 

-0.220 

- 0.014 

2 

C 

-0.063 

0.068 

-0.034 

-0.049 

3 

C 

0.055 

0.057 

-0.144 

- 0.001 

4 

C 

- 0.063 

- 0.087 

-0.165 

0.003 

5 

C 

0.207 

0.211 

- 0.036 

- 0.023 

6 

C 

- 0.070 

-0.060 

0.153 

0.017 

7 

C 

0.048 

0.078 

0.120 

- 0.001 

8 

0 

- 0.255 

-0.631 

-1.187 

0.002 

9 

C 

0.155 

-0.076 

0.233 

0.006 

10 

O 

- 0.353 

0.643 

-1.209 

0.029 

11 

0 

- 0.255 

-0.556 

1.157 

0.087 

12 

C 

0.005 

-0.039 

0.012 

-0.140 

13 

C 

-0.005 

-0.131 

0.141 

- 0.005 

14 

H 

0.013 

0.000 

0.000 

0.000 

15 

H 

0.004 

0.000 

0.000 

0.000 

16 

H 

-0.002 

0.000 

0.000 

0.000 

17 

H 

0.003 

0.000 

0.000 

0.000 

18 

H 

0.166 

0.000 

0.000 

0.000 

19 

H 

-0.016 

0.000 

0.000 

0.000 

20 

H 

-0.017 

0.000 

0.000 

0.000 

21 

H 

0.007 

0.000 

0.000 

0.000 

22 

H 

0.005 

0.000 

0.000 

0.000 

23 

H 

0.003 

0.000 

0.000 

0.000 

24 

H 

0.001 

0.000 

0.000 

0.000 

25 

H 

0.003 

0.000 

0.000 

0.000 


*e.u. = electron unit 
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Table 3 - Molecular charge distribution of p-n- octyloxybenzoic acid. 


Atom So. 

Symbol 

Charge 

(e.u.)* 

Atomic dipole componnts (Debye) 

X Y z 

I 

C 

0.393 

-0.220 

-0.221 

- 0.014 

i 

C 

- 0.063 

0.067 

- 0.034 

-0.049 

3 

c 

0.055 

0.057 

-0.144 

-0.001 

4 

c 

- 0.063 

-0.087 

-0.164 

0.005 

5 

c 

0.207 

0.209 

0.036 

- 0.025 

6 

c 

- 0.070 

-0.060 

0.152 

0.016 

7 

c 

0.048 

0.077 

0.120 

-0.001 

S 

0 

- 0.226 

- 0.637 

- 1.182 

- 0.059 

9 

0 

- 0.353 

0.643 

- 1.209 

0.029 

10 

0 

- 0.255 

- 0.556 

1.157 

0.087 

II 

c 

0.153 

- 0.086 

0.229 

0.022 

12 

c 

-0.001 

-0.044 

-0.127 

0.005 

13 

c 

0.022 

0.003 

0.115 

-0.017 

14 

c 

0.017 

- 0.025 

- 0.079 

0.021 

15 

c 

0.018 

0.030 

0.083 

- 0.022 

16 

c 

0.022 

-0.010 

- 0.080 

0.016 

17 

c 

0.025 

0.035 

0.084 

-0.010 

18 

c 

-0.009 

-0.155 

-0.045 

0.007 

19 

H 

0.166 

0.000 

0.000 

0.000 

20 

H 

0.002 

0.000 

0.000 

0.000 

21 

H 

0.013 

0.000 

0.000 

0.000 

22 

H 

0.004 

0.000 

0.000 

0.000 

23 

H 

-0.002 

0.000 

0.000 

0.000 

24 

H 

-0.017 

0.000 

0.000 

0.000 

25 

H 

-0.017 

0.000 

0.000 

0.000 

26 

H 

0.007 

0.000 

0.000 

0.000 

27 

H 

0.004 

0.000 

0,000 

0.000 

28 

H 

-0.007 

0.000 

0.000 

0.000 

29 

H 

-0.006 

0.000 

0.000 

0.000 

30 

H 

-0.007 

0.000 

0.000 

0.000 

31 

H 

-0.006 

0.000 

0.000 

0.000 

32 

H 

-0.008 

0.000 

0.000 

0.000 

33 

H 

-0.009 

0.000 

0.000 

0.000 

34 

H 

-0.009 

0.000 

0.000 

0.000 

35 

H 

-0.009 

0.000 

0.000 

0.000 

36 

H 

-0.009 

0.000 

0.000 

0.000 

37 

H 

-0.010 

0.000 

0.000 

0.000 

38 

7 A 

H 

0.001 

0.000 

0.000 

0.000 

39 

4 Pi 

H 

-0.003 

0.000 

0.000 

0.000 

4y 

H 

0.001 

0.000 

0.000 

0.000 

*e.u. - electron unit 
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TRANSLATION (A*) 



Fig. 3 - Variation of total stacking energy with respect to sliding (translation) dong itte 

(X). Energy is expressed in kcal/mole. (a) p-n-methoxybenzoic acid, (b) p-n- propoxybenzotc acid. W 


p-n-octyloxybenzoic acid. 
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Fig. 4 - Variation of stacking energy components with respect to sliding along the -tons L^^ybenzdc 
S Energy is expressed in kcal/mole. Notations are same as mentioned in Fig. 4(a). (a) p-n- methoxy 

acid (b) p-n-propoxybenzoic acid, (c) p-n- octyloxybenzoic acid. 
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Hg . s - 

to the long molecular axis(F). Energy ' acid. 

Fig. 4(a). (a) p-n-methoxybenzoic acid, (b) p-n-propoxybenzoic acia, W p J 






S.N. TIWARI and NITISH K. SANYAL 


66 


Table 4 - Total energy*, binding energy**, total dipole moment and its components for systems studied. 


Mole 

ccie Total energy 

Binding energy 


Dipole moment (Debye) 



(atomic unit) 

(atomic unit) 

MX 

PY 

Mz 

Total (p) 

P.MB 

- 1 18.13 

-9.46 

0.12 

-2.22 

0.03 

2.23 

PPB 

- 135.53 

- 11.97 

0.27 

-2.23 

0.09 

2.25 

FOB 

- 178.99 

- 18.21 

0.30 

- 2.27 

0.00 

2.29 

* Tot 

** Bi 
geom 

a! energy corresponds to the sum of atomic as well as electronic energies of all the constituents of the 

tale in the equilibrium geometry’. 

tiding energy of a molecule is the difference between the total energy of the equilibrium molecular 
eiry and the sum of the atomic energies of the constituent atoms. 

Table 

5 - Interaction energy components for stacking and in- 

expressed in kcal/mole. 

plane minimum energy configurations. Energy is 

Intern 

ction 

Stacking energy 



In-plane energy 


energy terms PMB 

PPB 

POB 

PMB 

PPB 

POB 


-0.02 

0.21 

0.64 

-0.15 

-0.22 

-0.19 

E 

** qmi 

-0.01 

0.31 

0.82 

-0.24 

-0.45 

-0.41 

E ■ • 

- 1.35 

-0.71 

1.03 

-0.01 

-0.11 

-0.10 


- 1.38 

-0.19 

2.49 

-0.40 

-0.78 

-0.69 

^pat 

-0.35 

-0.42 

-0.43 

-0.34 

-0.26 

-0.30 

E : 

*■* disp 

- 14.01 

- 14.13 

- 20.47 

-3.31 

-3.91 

-6.64 

F 

*-rep 

6.99 

5.96 

7.17 

1.28 

1.92 

3.01 

E 

**iQi 

-8.76 

-8.78 

-11.24 

-2.77 

-3.03 

-4.62 


energies are evaluated. All the possible stacked configurations have been examined by 
allowing rotations through 1 80° about the Y-axis, Y-axis or both axes. The results obtained, 
have been graphically shown in Fig. 2 (a,b,c). It is clear from this figure that variations in 
the total stacking energy with respect to interplanar separation, are similar for all the cases 
corresponding to the condition when both the stacked molecules are lying overlapped 
exactly one above the other and the complexes are stablized at 4.0 A'. 
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Fig. 6 - Variation of stacking energy with respect to rotational angle. Energy K 
are same as mentioned in Fig. 4(a). (a) p-n-methoxybenzoic acid, 
p-n-octyloxybenzoic acid. 
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T u r T tt: 4 - Total energy*, binding energy**, total dipole moment and its components for systems studied. 


MTec. 

,x To; dl energy 

Binding energy 


Dipole moment (Debye) 



j atomic unit) 

(atomic unit) 

Hx 

M-y 


Total (p) 

PMB 

- iIS.13 

-9.46 

0.12 

-2.22 

0.03 

2.23 

PPB 

- 135 53 

- 11.97 

0.27 

-2.23 

0.09 

2.25 

PV-B 

- 178.99 

- 1S.21 

0.30 

- 2.27 

0.00 

2.29 

* Tata! 
mL,ieeu! 

*-ik*uu 

gd-tlTLiitJ 

energy corresponds tc the sum of atomic as well as electronic energies of all the constituents of the 

!e in the equilibrium geometry, 

! T.g energy of a molecule is the difference between the total energy of the equilibrium molecular 
*\ and the sum of the atomic energies of the constituent atoms. 

Table 5 

- Interaction energy components for stacking and in- 
expressed in kcal/mole. 

plane minimum energy configurations. Energy is 


a:; 

Stacking energy 



In-plane energy 


energy t 

crui'- PMB 

PPB 

POB 

PMB 

PPB 

POB 

E . 

ri 

O 

© 

1 

0.21 

0.64 

-0.15 

-0.22 

-0.19 

r 

-0.01 

0.31 

0.82 

-0.24 

-0.45 

-0.41 


- 1 35 

-0.71 

1.03 

-0.01 

-0.11 

-0.10 

Ak 

-1.38 

-0.19 

2.49 

-0.40 

-0.78 

-0.69 

£ , 

-0.35 

-0.42 

-0.43 

-0.34 

-0.26 

-0.30 

£ 

- 14.U1 

- 14.13 

- 20.47 

- 3.31 

-3.91 

-6.64 


6.99 

5.96 

7.17 

1.28 

1.92 

3.01 


-8.76 

-8.78 

-11.24 

-2.77 

-3.03 

-4.62 


energies are evaluated. All the possible stacked configurations have been examined by 
allowing rotations through 180” about the X-axis, 7-axis or both axes. The results obtained, 
ha\e been graphically shown in Fig. 2 (a,b,c). It is clear from this figure that variations in 
me total stacking energy with respect to interplanar separation, are similar for all the cases 
corresponding to the condition when both the stacked molecules are lying overlapped 
exacth one above the other and the complexes are stablized at 4.0 A”. 
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Fig. 6 - Variation of stacking energy with respect to rotational angle. Energy is expressed ir. keaL mole. Notali' 
are same as mentioned in Fig. 4(a). (a) p-n-methoxybenzoic acid, (b) p-n- propo\\ benzoic acid, 
p-n-octyloxybenzoic acid. 
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However, the anisoptropic nature of pair interactions is reflected in the stacked 
complexes for PMB and PPB possessing maximum attraction at 3.5 and 4.0 A° respectively 
with a rotation of 180° about the X-axis while the minimum energy stacked complex for 
FOB corresponds to the overlapped structure at a separation of 4.0 A 0 (Fig. 2). Rotation 
about F-axis by 180° shows minimum attraction and that about both the axes gives 
intermediate attraction showing a minimum at 3.5 A° in all the cases. This implies that in 
dosed molecular packing, PMB and PPB are stacked such that in a molecular pair, one of 
the molecules is rotated by 180° about its own long molecular axis while no such rotation 
about any axis is needed for stablizing the molecular pair of POB molecules. It may be 
observed here that the effective range of the stacking interactions is of the order of 6.0 A 0 
- interplanar separation of the dimer (Fig. 2). 

Sliding : The translational freedom of one molecule over the other, is examined by 
sliding one molecule with respect to the other along both X and Y axes corresponding to 
the lower energy stacked molecular pair obtained from Fig. 2. For X-sliding (translation), 
four preselected rotations about the Z-axis viz . Z( 0°), Z(90°), Z(180°), Z(270°) were 
introduced with a view” to reach the extreme conditions and corresponding to each rotation, 
sliding energy is calculated in the range of ± 12.0 A° at an interval of 2.0 A°. The results 
have been plotted in Fig. 3 (a,b,c). As is evident from this figure, the sliding energy 
corresponding to Z{90°) and Z( 270°) has relatively low'er value. This is obvious because 
of the least molecular overlapping in these cases. A sharp and deep minimum is located 
in the case of PMB (Fig. 3a) w f hile relatively flat and wider regions are observed for the 
other tw o molecules. It is interesting to note that the variation of sliding energy corresponding 
to each rotation is similar in the case of PMB. The maximum energy difference lies within 
2 kcal/mole of dimer which indicates that at increased thermal agitations, the probabilities 
of resting in parallel and perpendicular stacking patterns become almost equivalent. This 
seems to be the reason why PMB directly goes to an isotropic melt state which represents 
a completely disorganised three dimensional structure. 

In [ be case of PPB, the two curves corresponding to Z(0°) and Z(180°) rotations show 
fiat regions. For Z(0 C ), the region starts from + 2.0 to - 4.0 A 0 , a span of 6.0 A° in which 
energy variations are minor (approximately 1 kcal/mole) while a similar flat region in the 
curve corresponding to Z(180°) is noticed between +4.0 to 1 1.0 A°, a span of 7.0 A\ These 
fiat regions indicate the existence of translational freedom of the molecules along their 
long molecular axes. The energy for Z(180°) is slightly lower than that of Z(0°) and as 
such, has been undertaken for further study. In the case of POB, flat regions are obtained 
for both Z(0 ) and Z(180 ), but the energy corresponding to Z(180°) is higher as compared 
to Z(0 ) and, hence, the minimum located by Z(0°) curve, is preferred for further calculations. 
Rapid increase in energy, observed between 3.0 and 8.0A° appears due to the violation of 
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Fig. 7 - Minimum energy stacked complexes for (a) PMB, (b) PPB and (c) POB molecules with interplanar 
separation (a) 3.3A\ (b) 3.6A’ and (c) 3.7 A 0 respectively. Energy details are given in Table 5. 



Fig. 8- Minimum energy in-plane configurations for (a) PMB, (b) PPB and (c) POB molecules with 
intermolecular separation (a) 6.1 A', (b) 6.5A* and (c) 7.0 A" respectively. Energy details are given m 
Table 5. 
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van der Waal contact distances for C(H)— C(H) between the methyl groups. It may be 
noted that the stacked pair energy is minimised at 4.0 A° (Fig. 2c) while C(H) — C(H) 
contact distance is 4.20 A 0 . Hence, by slightly increasing the interplanar separation, one 
may sufficiently reduce such repulsive forces. 

The variation of different interaction energy components corresponding to the rotational 
minima obtained from Fig. 3 (a,b,c), has been shown in Fig. 4 (a,b,c). As is evident from 
this figure, polarization component is weak while electrostatic energy curve assumes both, 
attractive and repulsive nature but is insignificant. The only dominant term is due to ‘short 
range’ dispersion and repulsion components. The repulsion energy term has not been plotted 
explicitly as it can be easily obtained from Kitaigorodskii (dispersion + repulsion) and 
dispersion energy curves. The repulsion region, noticed for POB, is also supported by 
Kitaigorodskii curve implying that though dispersion energy is highly attractive but the 
sudden increase in the repulsive energy results in the repulsive nature of the total energy. 
The dominant role of dispersion forces (induced dipole-induced dipole) in stabilizing the 
packing of mesogens, has been also recognized by several other workers 12,44,45 . 

The variation of intermolecular interaction energy components for sliding along T-axis, 
has been shown in Fig. 5 (a,b,c) respectively for PMB,PPB and POB molecules. It is 
obvious from this figure that only for translations within the range of ± 2.0 A\ energy 
largely changes in case of PMB. In other words, a change of approximately 4.5 kilocal/mole 
of dimer is observed for PMB while PPB and POB exhibit relatively smaller changes in 
the energy (approximately 2 kcal/mole). However, this change in energy is higher as 
compared to sliding along X-axis as shown in Fig. 3 (a,b,c) implying that there is a lesser 
possibility oi sliding of the stacked molecules along an axis perpendicular to the long 
molecular axis in the molecular plane. Gross similarity exist in the total and Kitaigorodskii 
(dispersion + repulsion) energy curves. 

Rotation : The variation of interaction energy components with respect to rotation 
corresponding to the minimum energy points obtained from Fig. 5 (a,b,c) has been plotted 
in Fig. 6 (a,b,cj respectively for PMB, PPB and POB molecules. Both, clockwise and 
anticlockwise rotations in the range of ± 6Q\ has been given with an aim to examine the 
angular dependence of stacking energy. It is interesting to note from Fig. 6(c) that maximum 
selectivity is obtained for POB while PMB shows maximum orientational freedom. The 
orientational freedom of PPB is intermediate with minor changes in energy (approximately 
1 kcal/mole) within range of ± 20°. Also the minimum for POB corresponds to -10.0 
kcal/mole which is largest in magnitude as compared to PMB or PPB molecules. This 
indicates that the tendency of forming aligned structure is maximum in the case of POB 
which supports the smectic character of system at lower temperatures in the mesomorphic 
range. 
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The minimum energy configurations were subjected to further refined calculations to 
achieve the final lowest energy configurations. The minimum energy stacked complexes 
thus obtained, have been shown in Fig. 7 (a,b,c) while the energy distribution has been 
listed in Table 5. It is observed that two molecules of PMB are stabilized at 3.30 A° such 
that one of them has been rotated by 180° both about X and Z-axes. It is interesting to note 
that PMB and PPB both have the same order of energy values (nearly 9.0 kcal/mole) but 
differ in their respective stacking patterns. Though the stacking pattern of POB molecules 
is similar to PMB, there is an increase in interplanar separation (3.70 A°) with a lowering 
in the stacking energy (-11.24 kcal/mole for POB) which is probably due to the long 
alkoxychain attached to POB. The minimum energy configurations obtained from in-plane 
interactions, are shown in Fig. 8(a,b,c). It is significant to note that all the homologues 
under study exhibit similar in-plane minimum energy configurations. The in-plane energy 
for PMB and PPB are nearly the same - 3.0 kcl/mole while for POB the energy corresponds 
to - 4.62 kcal/mole. This implies that by increasing the alkoxy chain, in-plane attractions 
are largely increased as compared to the stacking which is in accordance with the 
experimental observations and likely to favour the smectic behaviour of the compound 1 . 

Conclusion 

Computational results based on the consideration of actual molecular shape and 
composition indicate that the packing energy of the members of p-n-alkoxybenzoic acid 
series is in the following order : 

POB > PPB > PMB 

This shows that the stability of molecular packing increases with the addition of 
methylene group in the alkoxychain. The dominant forces between the mesogenic molecules 
are essentially of dispersion (induced dipole-induced-dipole) type which is in agreement 
with the basic assumption of the molecular field theory. Further, the analysis of energy 
variation with respect to orientation and sliding is helpful in discriminating mesogenic and 
non-mesogenic molecules. These studies may also be useful in characterizing smectic and 
nematic properties of molecules. 
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Abstract 

Dielectric relaxation study of pyridine-alcohol viz. methanol, ethanol, propan-l-oi and butan-l-ol 
mixtures have been carried out for 1 1 different concentrations at 25 °C. Time Domain Reflectometiy 
(TDR) in reflection mode has been used to measure reflection coefficient in frequency range of 10 
MHz to 10 GHz. The bilinear calibradon method has been used to obtain dielectric parameters viz. 
static dielectric constant and relaxation time. The concept of effective volume is introduced and the 
values of effective volume occupied by the pyridine in presence of alcohols are determined by using 
four models, based on the excess permittivity, excess relaxation time, the Bruggeman equation and 
the Kirkwood equation. The values of a parameter a which represents change in volume of constituent 
1 due to constituent 2, in the mixture of constituents 1 and 2, are also determined using different 
models. The volume shrinks for pyridine in presence of methanol and ethanol, whereas this expands 
m presence of butan-l-ol and propan- l-ol. 

(Keywords : time domain reflectometry/dielectric parameters/alcohols/pyridine) 

Introduction 

The dielectric relaxation study of solute-solvent mixtures at microwave range gives 
information about the molecular polarization in the system. This is related to formation of 
multimers in dipoles and their rotations. 

The present paper reports the dielectric study related to pyridine and alcohol mixtures. 
Methanol, Ethanol, Propan-l-ol and Butan-l-ol are used as alcohols. The study was done 
at 25 C. The Time Domain Reflectometry in reflection mode has been used for the study 1 * 5 . 

Materials and Method 

Pyridine (LR grade). Methanol, Ethanol (Spectroscopic grade), Propan-l-ol and 
Butan-l-ol (AR grade) were obtained commercially and used without further purification. 
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The solution were prepared at different volume percentage of pyridine (t^) in alcohols in 

step of 107^, at 25*C. 

Apparatus used in this study was the Tektronix 7854 sampling oscillascope with 7S12 
TDR unit. In this apparatus, a fast rising step voltage pulse of 200 mV amplitude and 25 
ps rise time with repetition frequency of 1 kHz is generated by a tunnel diode and is 
propagated through a coaxial transmission line. The sample is placed at the end of the 
coaxial transmission line in a Standard Military Application (SMA) coaxial cell. The SMA 
cell used for this work had 3.5 mm outer diameter and 1.35 mm effective pin length. The 
step pulse generated by tunnel diode and the pulse which is reflected from the sample cell 
were sampled by sampling oscilloscope in the time window of 5 ns. The reflected pulse 
without sample i?, (?) and with sample R x (?) were digitised into 1024 points in the 
oscilloscope memory and then transferred to PC/XT through GPIB. 

A typical example of the reflected pulses without sample and with sample is shown 

in fig. 1. 

The temperature controller system with water bath and a thermostat has been used to 
maintain the constant temperature within the accuracy limit of ±1°C. The sample cell is 
surrounded by a heat insulating container through which the water of constant temperature 
using a temperature controller system is circulated. The temperature at the cell is checked 
using the electronic thermometer. 


Data Analysis 

The time dependent data were processed to obtain complex reflection coefficient spectra 
p*(co) over the frequency range from 10 MHz to 10 GHz using Fourier transformation 6,7 

as 


P* (to) = (c/turf) [p (o >)lq (co)] (1) 

where p(w) and q( to) are Fourier transforms of (i? f (/) - R x (t)) and (R^t) + R x (t)) 
respectively, c is the velocity of light, co is angular frequency, d is effective pin length. 

The complex permittivity spectra e*(to) were obtained from reflection coefficient 
spectra p*(co) by using bilinear calibration method 8 . Alcohols and pyridine are chosen to 
be as calibrating liquids. 


The experimental values of t* are fitted with the Debye expression 9 ' 11 
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with e pt and r fn as fitting parameters. A non-linear least-squares fit method 12 was used to 

determine the values of dielectric parameters. 

The value of is taken to be fixed and computed as 


£oc, m = *1 + (1 ~ *i) e-o 2 

where e^r and e x2 are permittivities at infinite frequencies for constituents 1 and 2, 
respective!}', and x : is the volume fraction of the constituent 1 in the mixture. They are 

estimated to be as 3.5 and 4, respectively. 

Results and Discussion 

The density and e m values of pure liquids along with corresponding literature values 
are .given in Table 1, which shows a good agreement. The static dielectric constant and 
relaxation time obtained by the least squares fit method are listed in Tables (2 and 3). 

The solute-solvent interaction can be understood by considering non-linear effect in 
the values of dielectric parameters in the mixtures. The recent studies 13,14 have shown that 
the static dielectric constant (polarization) of associated liquids depend linearly on volume 
fraction. On the other hand, mixtures of associated liquids have shown relaxation times 
always larger than average of the individual rate of the pure components. The conventional 
approach in the case of ideal solutions of associated liquids with no significant structural 
changes with respect to pure components is that a linear behaviour of the logarithm of the 
relaxation rate should be obtained i.e. 

m 

ln(r(x lm )) = x lm In (Tj) + (i -jfiJ ln(T,) (3) 

where x lm is mole fraction of the component 1. m, 1 and 2 represent mixture, component 

! and 2 respectively. The suffixes 1 and 2 represent pyridine and alcohols, respectively, in 

the present paper. 

It is useful to define the effective volume fraction as follows ; 
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eff ('U) (®) 

4 = (4) 

and the effective mole fraction as 

eff _ In (x w (£ij) ~ In (x,„ (0)) 

Xlm ln(t m (l)) - ln(x ra (0)) (5) 

where e m (x } ) is the dielectric constant of mixture with the volume concentration of 
component 1 as x { in the mixture 1. 

The value of the x e f v may be calculated from xf m as follows : 



where/ = (p 2 M { / p { M 2 ) ; p- and M t represent density and molecular weight respectively, 
for constituent i (/ = 1, 2) . 

Table 1 - Comparison of data for the pure liquids used with the literature values at 25°C. 


Liquid 

This work 

<H) 

Lit. 

This work 

-3 

p/gm.cm 

Lit. 

Methanol 

32.05 

32.63 

0.78695 

0.78698 

Ethanol 

24.08 

24.30 

0.78509 

0.78506 

Propan- 1-01 

20.21 

20.10 

0.78124 

0.7812 

Butan-1-01 

17.22 

17.10 

0.80596 

0.8059 

Pyridine 

12.02 

12.30 

0.98113 

0.9811 


The plots of xf v vs x x are given in Fig. (2 and 3). The ideal behaviour is shown by the 
straight line. The effective values are less than the ideal behaviour for smaller alcohols. As 
size of the alcohol molecule increases, the effective volume increases. This indicates that 
there is a tendency to occupy more volume by pyridine molecules in presence of larger 
alcohols whereas it occupies smaller volume in presence of alcohol with smaller size. This 
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:;c dxlsctric parameters (£ m ) for pyridine-alcohol system at 25“C. 


Methanol 

3105 ( 11 ) 

28 . 87 ( 14 ) 

27 . 45 ( 8 ) 


25 . 27 ( 7 ) 


21 . 31 ( 9 ) 

19 . 78 ( 5 ) 

17 . 41 ( 9 ) 

15 . 14 ( 5 ) 

12 . 02 ( 2 ) 


24 . 08 ( 5 ) 

23 . 17 ( 10 ) 

21 . 73 ( 4 ) 

21 . 39 ( 4 ) 

20 . 11 ( 4 ) 

18 . 36 ( 3 ) 

17 . 26 ( 3 ) 

16 . 52 ( 4 ) 

14 . 84 ( 2 ) 

13 . 44 ( 4 ) 

12 . 02 ( 2 ) 


Propan- l-ol 
20 . 21 ( 1 ) 
17 . 22 ( 8 ) 
16 . 72 ( 16 ) 
15 . 68 ( 10 ) 
15 . 11 ( 5 ) 
14 . 45 ( 1 ) 
14 . 21 ( 4 ) 
13 , 86 ( 15 ) 
13 . 46 ( 21 ) 
12 . 92 ( 21 ) 
12 . 02 ( 2 ) 


Butan-l-ol 

17 . 22 ( 2 ) 

15 . 73 ( 11 ) 

15 . 12 ( 10 ) 

14 . 53 ( 5 ) 

14 . 23 ( 3 ) 

13 . 93 ( 8 ) 

13 . 49 ( 3 ) 

12 . 89 ( 4 ) 

12 . 71 ( 4 ) 

12 . 38 ( 2 ) 

12 . 02 ( 2 ) 


means that the volume occupied by pyridine exerts larger effective repulsive forces to 
volume occupied by the alcohols, as the size of the alcohol increases 

In the paper 4 , the Bruggeman formula was modified to exDlain th* , 

experimental data related to binary polar liquids as follows : ‘° ar of the 


e mC*i) ~ UUl fM2) 


1><2) ~ ' £ m (1) J I £ m (1) 


l-jf 


x w = [a-(u- l)xj]x, 


The above equation satisfies the boundary conditions at x - o „ j _ 

^ r ^ __ i , . . i "" u ancl — i. I he value 

ot a - I corresponds to no interaction between constituent 1 anH -> „ ^ , 

reduction of the effective volume occupied by the constituent l h. ’ corresponds to 

» > > fronds to incase in *e effect volume conalZlT' 1 * T* 
constituent 2, onstituent 1, due to the 
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Table 3 - Relaxation time (t m ) for pyridine-alcohol system at 25”C. 


x l 

Methanol 

Ethanol 

Propan- l-oi 

Butan-l-ol 

00 

52.37 

145.84(80) 

278 86*11) 

512.04(17) 

0.1 

48.07(45) 

102.46(10) 

148.43(14) 

212.35(52) 

0.2 

44.82(27) 

94.28(43) 

110.42(18) 

141.72(31) 

0.3 

42.36(17) 

73.18(31) 

78.90(11) 

103.98(10) 

0.4 

34.87(25) 

56.54(19) 

59.82(82) 

85.13(21) 

0.5 

30.63(3) 

45.37(39) 

41.12(4) 

61.60(13) 

0.6 

28.20(32) 

33.62(21) 

35.77(25) 

49.42(33) 

0.7 

25.02(24) 

27.77(19) 

28 43(91) 

36.76(38) 

0.8 

19.45(32) 

24.22(15. 

22.17(47) 

27.21(36) 

0.9 

16.75(24) 

18.82(82) 

17.75(32) 

19.49(28) 

1.0 

10.11(29) 

10.11(29) 

10.11(29) 

10.11(29) 

Table 4 - 

Values of a for pyridine-alcohol systems. 



‘ a' 

Methanol 

Ethanol 

Propan- 1 -01 

Butan'l-ol 

a eQ 

0.770 

0.780 

1.719 

1.791 

a T 

- 0.450 

0.371 

1.309 

1.107 

a B 

0.494 

0.582 

1.606 

1.626 

a K 

0.901 

0.914 

1.652 

1.756 

The 

Kirkwood correlation 

factor g 15 is also a 

parameter for getting structural 


information. In the mixture, the effective volume xffi can be obtained by using the following 


Kirkwood equation for mixtures, 


4YIN 

9kT 


ViPigj 

M, 




Mv 


Af, 



(£>n ^om) 

£ ,'p + 2) 2 


V 


(9) 
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where 1,2) are dipole moments of constituents 1, 2; g : (/ 1,2) are correlation 

factors of constituents 1, 2; k is the Boltzmann factor and T is temperature of the mixture. 

The correlation factors of pure liquids 1 and 2 are computed by using the Kirkwood 
equation for pure liquids. By using Eqn. (8) and (9), it is possible to determine the 

corresponding \alues of a. 

t a z , Qg and are values oi ci as determined by values of Ay v from Eqn. (4), (5 ), 
(7) and (9), respectively, and using Eqn. (6) and (8). These values are listed in Table 6. It 
is interesting to observe that as size of alcohol increases, the value of a also increases. This 
indicates that the effective volume of pyridine gets increased due to the larger size of the 

alcohols. 

The valuo of a obtained from the relaxation data are much smaller than corresponding 
values of a obtained from polarization volume. This is expected, as relaxation time involves 
dynamical behaviour of polarization, so more sensitive to long range intermolecular 
interactions. 


Conclusion 

The Onsagar dielectric of pure liquids assumes a cavity inside the liquid. The radius 
of the cavity is computed by using the density. We have computed the values of effective 
volumes in the mixture by using four different models based on the excess permittivity, 
excess logarithmic of relaxation time, the modified Bruggemann equation and the Kirkwood 
equation, respectively. The parameter a which represents significance of change in volume, 
is determined. This provides a unified approach to understand the interaction between two 
constituents in the liquid mixtures, using different models. 

Acknowledgement 

Financial support from Department of Science and Technology, New Delhi, is gratefully 

acknowledged. 

References 

1. Kumbharkhane, A.C., Helambe, S.N. & Mehrotra, S.C. (1993) J . Chem. Phys, 23 : 2205. 

2. Puranik, S.M., Kumbharkhane, A.C. & Mehrotra, S.C (1992) J. Chem, Soc . Faraday Trans. 88 : 433. 

3. Puranik, S.M., Kumbharkhane, A.C. & Mehrotra, S.C (1991) J. Chem. Soc. Faraday Trans. 87 : 1569. 

4 Kumbharkhane, A.C., Puranik, S.M. & Mehrotra, S.C. (1991) J. Sol. Chem. 20 : 12. 



DIELECTRIC STUDY OF PYRIDINE-ALCOHOL BINARY LIQUIDS 83 

5. Pui-mk, S.M., Kumbharkhane, A.C. & Mehrotra, S.C. (1991) 7. Micro. Pirn. EM energy 26 : i%. 

6. Shannon C.E. (1949) Proc. IRE 37 : 10. 

7 Sarnulon, H. A. (1951) Proc. IRE 39 : 175 

g. Cole, R H , Berbarian, J G.,Mashimmo,S.,Chr>ssiko$.G.. Burns, A.&Tomban E. (1989 ,7 -W 
66:793. 

9. Havnlliak, S. & Negami, S. (1996) 7: Po/ym. Set. Po/ym. Svmp. 99. 

10 Cole, K S. & Cole, R.H (1941) 7. Chem. Phys. 9 : 341. 

I i Davidson, D.W. & Cole, R.H. (1950) J. Chem . Phys. 18 : 1484 

12 Bevington, P.R. (1969) Data reduction and error analysis for the physical sciences. Me G raw Hill. New* 

York. 

13 Mashimo, S., Umehara, T. & Redlm, H. (1991)7. Chem. Phys. 95 : 62~’5. 

14. Bao, J., Swicord, M.L & Davies, C.C. (1995) 7. Chem. Phys. 104 : 4441 

15. Frohlich, H. (1949) Theory of Dielectrics, Oxford University Press, London. 



1 DITORJAL BOARD 


Chief Editor 
Prof. H.C. Khare 


Former. Professor of Mathematics, University ot Allahabad. 
The National Academy of Sciences. India, 

5 f Lajpatrai Road, Allahabad - 21 1 002 
Fax ; 091-532-641 ! S3 
E-Mail i ruisi-c .ide.vsnl.iicl.in 


Prcf. R.P. Agrawal 
Former, Vice-Chancellor, 

Rajasthan & Lucknow Universities, 
81/201, Nirala Nagar, 

Lucknow - 226 020 
(Mathematics) 


Prcf. Afak K. Gupta 

Head, Department of Earth and 

Planetary Sciences, 

University of Allahabad, 

Allahabad-211 002 
Fax: 091-532-641340 
E-Mail : ncemp@nde.vsnl.net in 
(Earth Sciences) 

Prof. E.L. Khandeiwai 
Emeritus Scientist (CSSR), 

Defence Materials and Stores Research 
and Development Establishment, 
DMSRDE Post Office. 

G.T. Road, 

Kanpur - 208 013 

Fax:091-512-450404 

(Chemistry) 


2. Prof. Suresh Chandra 
Department of Physics, 

Banaras Hindu University, 

Varanasi - 221 005 

Fax '091-542-317074 

E-Man' . sJumdraui banaras cmct.in 

(Physics) 

4. Prof. S.V Kessar 
Senior Scientist, 

Department of Chemistry, 

Panjab University, 

Chandigarh - 160 014 

E-Mail : sv kessaffo?paniabum%.c h<.i n e . m 

(Chemistry) 


6. Prof. U.C. Mohanly 
Professor & Head, 

Centre for Atmospheric Science 
Indian institute of Technology, 
Hauz Khas, 

New Delhi - 110 016 
Fax:091-11-6862037 
E-Mail : mohantvAi cas.iitd.crnct.m 
(Atmospheric Sciences) 


Prof.. K.S, Valdiya 
Bhatnagar Research Professor, 

Jawaharial Nehru Centre for 
Advanced Scientific Research, 

Jakkur P.0, 

Bangalore - 560 064 
Fax * 091-80-6462766; 

E-Mail : udaviiFincasi.ac.in 
(Geology) 

Managing Editor 
Prof. S.L. Srivastava 

Formerly Professor & Head, Department of Physics, University of Allahabad; 

• The National Academy of Sciences, India, 

5, Laipatrai Road, Allahabad - 21 1 002 
Fax:091-532-641183 
E-Mail : nasrO'-nde.vsnl.neUn 



EDITORIAL ADVISORY BOARD 


1 . prof. Edwin D. Becker 
Scientist Emeritus, 

National Institutes of Health, 

5 Center Drive, 

Bethesda, 

Maryland 20892-0520, U.S.A. 

E-Mail :tbecker@mh.gov 

(Spectroscopy/Nuclear Magnetic Resonance) 

3. Prof. S. Chandrasekhar 
Director, 

Centre for Liquid Crystal Research, 

P.B. No.1329, 

Bangalore - 560 013 
Fax : 091-80-8382044 
E-Mail:uc!er@giasbg01.vsnl.net.in 
(Liquid Crystals/Condensed Matter) 

5. Prof. M.G.K. Menon 

President, ISI & Past President ICSU, 

Chairman, Board of Governors, HT (Bombay), 
C-63, Tarang Apts., 

19-I.P. Ext., Mother Dairy Road, 

Patparganj, 

Delhi- 110 092 

Fax:091-11-6959456 

E-Mail : mgkmenon@.ren02.nic.in 

(Physics) 


7. Prof. C.Kumar N. Patel 
Chairman & CEO, 

Pranalytica Inc., 

1101 Colorado Avenue, 

Santa Monica, CA90401 - 3009, U.S.A. , 

Fax:310-458-0171 

(Physics) 


9. Dr. P. Rama Rao 
Vice-Chancellor, 

University of Hyderabad, 

P.0. Central University, 

Hyderabad - 500 046 

(Physical and Mechanical Metallurgy) 


11. Prof. Govind Swarup 
Homi Bhabha Senior Fellow, 

National Centre for Radio Astrophysics, 
Tata Institute of Fundamental Research, 
Pune University Campus, 

P. B. No. 3, 

Ganeshkhind, 

Pune -411 007 

Fax : 091-20-5657257, 5655149 

E-Mail : eswarup@ncra.tifr.res. in 

(Radio Astronomy) 


2. Prof. Sir Harmann Bondi 
60 Mill Lane, Impington, 
Cambridge CB4 9XN, U.K. 
Tele: 01223 565180 
(Mathematical Astronomy) 


4. Prof. S.K. Joshi 

Sarabhai Research Professor, JNCASR, 

National Physical Laboratory, 

Dr. K.S. Krishnan Marg, 

New Delhi -110 012 
Fax:091-11-5852678 
E-Mail : skjoshi@csnpI.ren.nic.in 
(Solid State Physics) 

6. Prof. A.P. Mitra 

Formerly Secretary, DSIR and 
Director-General, CSIR, 

Hon. Scientist of Eminence, 

National Physical Laboratory, 

Dr. K.S. Krishnan Marg, 

New Delhi -110 012 
Fax:091-11-5852678, 5764189 
E-Mail : 3pmitra@doe.emet.1n : apmitra@ndf.vsnl.netm 
(Ionospheric Physics/Radio Communication/ 
Space Physics/Space Science) 

8. Dr. B.L.S. Prakasa Rao 
Distinguished Scientist, 

Indian Statistical Institute, 

7, S.J.S. Sansanwal Marg, 

New Delhi - 110 016 
Fax:091-11-6856779 
E-Mail : bisp@isid.ac.in 
(Mathematical Statistics) 

10. Prof. M.M. Sharma 

Kothari Research Professor (Hon.), JNCASR, 

Formerly Director & Professor of Chemical 
Engineering, UDCT, 

502, Saurabh, Plot No. 39, 

Kunder Marg, Swastik Part, Chembur, 

Mumbai - 400 071 

E-Mail : mmsharrna@b nm-3.vsnl.neUs 

12. Prof. H.C. Khare (Chief Editor) 

Former, Professor of Mathematics, 

University of Allahabad, 

The National Academy of Sciences, India, 

5, Lajpatrai Road, 

Allahabad -211 002 

Fax:091-532-641183 

E-Mail : nasi@nde.vsnl.netm 

{Applied Mathematics/Theoretical Physics] 



